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Abstract 



O , 
(D . 

Q ■ A complete and critical survey is presented of all half-life, decay-energy and branching-ratio 

*^ . measurements related to 20 superallowed 0"*" — > 0"*" decays; no measurements are ignored, though 

I \ some are rejected for cause and others updated. A new calculation of the statistical rate function / 

' is described and experimental ft values determined. The associated theoretical corrections needed 

in 

• to convert these results into J^t values are discussed, and careful attention is paid to the origin 

<N ; 

^—1 , and magnitude of their uncertainties. As an exacting confirmation of the conserved vector current 

' hypothesis, the J-t values are seen to be constant to 3 parts in 10^. These data are also used to 

• set a new limit on any possible scalar interaction: Cs/Cy = —(0.00005 it 0.00130). The average 
^ . J^t value obtained from the survey, when combined with the muon liftime, yields the up-down 

, quark-mixing element of the Cabibbo-Kobayashi-Maskawa matrix, Vud = 0.9738 it 0.0004; and the 

' unitarity test on the top row of the matrix becomes |T4rfp -|- |T4sp -|- |T4bP = 0.9966 ± 0.0014 using 

■ the Particle Data Group's currently recommended values for Vus and Vub- We also express this 

result in terms of the possible existence of right-hand currents. Finally, we discuss the priorities 
for future theoretical and experimental work with the goal of making the CKM unitarity test more 
definitive. 
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I. INTRODUCTION 



Precise measurements of the beta decay between nuclear analog states of spin, J'^ = 0"*", 
and isospin, T = 1, provide demanding and fundamental tests of the properties of the 
electroweak interaction. Collectively, these transitions can sensitively probe the conservation 
of the vector weak current, set tight limits on the presence of scalar or right-hand currents and 
contribute to the most demanding available test of the unitarity of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix, a fundamental tenet of the electroweak standard model. 

Eight transitions, ^^O, ^eAP, 34q^ ssj^m^ 42gp^ m^j ^ so^n and ^^Co are particularly 
amenable to experiment and, because of their significance to physics, have consequently 
received a good deal of attention over the past few decades. In each of these cases, the 
experimental /t- value is known to better than 0.1%. In the 1990s, was added to this 
list; its ft value is known to a precision of 0.15%. More recently, three more cases have been 
added: ^^Mg, ^^Ar and ^^Rb, with /t-value standard deviations ranging from from 0.24% to 
0.40%. In the near future these uncertainties will undoubtedly be reduced and an additional 
eight cases could well be added to the list. Though improvements are still possible, with 
current data we can test the conserved vector current hypothesis at the level of 3 parts in 
10^ and the three-generation Standard Model at the level of its quantum corrections. 

Over the past decade, it has become increasingly clear that the CKM unitarity test 
made possible by these measurements does not, in fact, quite agree with standard- model 
expectations. The test involves the top row of the CKM matrix and requires that the sum 
of squares of the three experimentally-determined elements, |Ktdp + |V"„sp + iKtfePj should 
equal 1. With results from superallowed /3-decay providing the input for Vud, and values for 
Vus and Vub taken from the Particle Data Group reviews, the sum falls short of unity by 0.3 
%, more than twice the quoted standard deviation 1] - a provocative but hardly definitive 
disagreement. Nevertheless, it has stimulated experimental activity not only on the nuclear 
decays used to determine Vud but also on the i^es. branching ratio used for Vus- Strikingly, 
a new measurement of the K'^^^ branching ratio [2j has thrown the accepted value of Vus 
into doubt. Although the new branching-ratio result disagrees significantly with previous 
measurements, it would, if taken by itself, lead to a larger value for Vus and thus bring 
the CKM top-row sum into agreement with unity. At this time, the value of Vus remains 
controversial and there are a number of kaon-decay experiments currently underway, which 
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should lead to a settled outcome within a very few years. 

With all this activity in progress, and the likelihood that a new and reliable value of 
Vus will soon be forthcoming, this is an opportune time to produce a complete new survey 
of the nuclear data used to establish Vud- This way, we will be able to view the value of 
Vud with renewed confidence in anticipation of a revised result for V^s- {Vut is very small 
and contributes a negligible .001% to the unitarity sum.) We have published four previous 
surveys, refs. Bflflfl: the most recent appeared fifteen years ago and included only eight 
superallowed transitions. In addition to bringing the results for these cases up to date, we 
are now incorporating data on twelve more transitions and have continued the practice we 
began in 1984 of updating all original data to take account of the most modern calibration 
standards. We have also made completely new calculations of the statistical rate function, 
/, and employed the most complete radiative and isospin-symmetry-breaking corrections in 
dealing with the /t- values in the context of fundamental weak-interaction tests. 

Superallowed Fermi beta decay between 0+ states depends uniquely on the vector part 
of the hadronic weak interaction. When it occurs between isospin T = 1 analog states, the 
conserved vector current (CVC) hypothesis indicates that the ft values should be the same 
irrespective of the nucleus, viz. 

ft = = constant, (1) 

where K/{hcf = 27t^h\n2/ {niec'^f = (8120.271 ± 0.012) x 10"^° GeV~^s, Gy is the vector 
coupling constant for semi-leptonic weak interactions, and Mp is the Fermi matrix element, 
which for T = 1 states has the value Mp = \f2. The CVC hypothesis asserts that the 
vector coupling constant, Gv, is a true constant and not renormalised to another value in 
the nuclear medium. A demonstration with the data assembled here that the ft values are 
indeed constant would provide a stringent test of the CVC hypothesis. 

Unfortunately, Eq. has to be amended slightly. Firstly, there are radiative corrections 
because, for example, the emitted electron may emit a bremsstrahlung photon, which goes 
undetected in the experiment. Secondly, isospin is not an exact symmetry in nuclei so the 
nuclear matrix element, Mp is slightly reduced from its ideal value, leading us to write: 
iMiT'P = 2(1 — 5(7). Thus, we define a "corrected" ft value as 

Tt = ft{\ + bR){\ - be) = 20^1 + ^ (2) 
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where 5c is the isospin-symmetry-breaking correction, 5r is the transition-dependent part 
of the radiative correction, and is the transition-independent part. Fortunately these 
corrections are all of order 1% but, even so, to maintain an accuracy criterion of 0.1% they 
must be calculated with an accuracy of 10% of their central value. This is a demanding 
request, especially for the nuclear-structure-dependent corrections. 

To separate out those terms that are dependent on nuclear structure from those that are 
not, we split the transition-dependent radiative correction into two terms. 



of which the first, 5^, is a function only of the electron's energy and the charge of the 
daughter nucleus Z] it therefore depends on the particular nuclear decay, but is independent 
of nuclear structure. The second term, 5ms^ like 5c, depends in its evaluation on the details 
of nuclear structure. To emphasize the different sensitivities of the correction terms, we 
rewrite the expression for as 



where the first correction in brackets is independent of nuclear structure, while the second 
incorporates the structure-dependent terms. 

Our procedure in this paper will be to examine all experimental data related to 20 super- 
allowed transitions, comprising those that have been well studied, together with others that 
have only recently become accessible to precision measurement. The methods used and the 
data accepted are presented in Sect.lHl The calculations and corrections required to extract 
final jFt- values from these data are described and applied in Sect. in the same section, 
we use the resulting jFt- values to test CVC Finally, in Sect. IIVI we explore the impact of 
these results on a number of weak-interaction issues: CKM unitarity as well as the possible 
existence of scalar or right-hand currents. 

II. EXPERIMENTAL DATA 

The /t-value that characterizes any /^-transition depends on three measured quantities: 
the total transition energy, Qeci the half-life, ti/2, of the parent state and the branching 
ratio, i?, for the particular transition of interest. The Qec- value is required to determine 



5r = 5'ji + 5ns, 



(3) 



J^t = ft{l + 5'j,){l + 5Ns-5c), 



(4) 
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the statistical rate function, /, while the half-life and branching ratio combine to yield the 
partial half-life, t. In tables ITl lVIII we present the measured values of these three quantities 
and supporting information for a total of twenty superallowed transitions, incorporating the 
eight cases we dealt with in our last complete survey [6] , but now including four more cases 
that have been measured more recently with comparable precision, and a further eight that 
are likely to become accessible to precision measurements within the next few years. 



A. Evaluation principles 

In our treatment of the data, we considered all measurements formally published before 
November 2004 and those we knew to be in an advanced state of preparation for publication 
by that date. We scrutinized all the original experimental reports in detail. Where necessary 
and possible, we used the information provided there to correct the results for calibration 
data that have improved since the measurement was made. If corrections were evidently 
required but insufficient information was provided to make them, the results were rejected. 
Of the surviving results, only those with (updated) uncertainties that are within a factor of 
ten of the most precise measurement for each quantity were retained for averaging in the 
tables. Each datum appearing in the tables is attributed to its original journal reference via 
an alphanumeric code comprising the initial two letters of the first author's name and the 
two last digits of the publication date. These alphanumeric codes are correlated with the 
actual reference numbers in Table lyTTTl 

The statistical procedures we have followed in analyzing the tabulated data are based 
on those used by the Particle Data Group in their periodic reviews of particle properties, 
e.g. ref and adopted by us in earlier surveys J, 0| of superallowed 0+ —>■ 0+ beta decay. 
In the tables and throughout this work, "error bars" and "uncertainties" always refer to 
plus-and-minus one standard deviation (68% confidence level). For a set of uncoupled 
measurements, Xi ± 6xi, of a particular quantity, a gaussian distribution is assumed, the 
weighted average being calculated according to: 

x±Sx=^^±{E^wr'^\ (5) 

where 

Wi = 1/ {5xi) 
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and the sums extend over all measurements. For each average, the is also calculated 
and a scale factor, S, determined: 



This factor is then used to establish the quoted uncertainty. If 5 < 1, the value of 6x from 
Eq. (jni) is left unchanged. If S" > 1 and the input 6xi are all about the same size, then 
we increase 6x by the factor S, which is equivalent to assuming that all the experimental 
errors were underestimated by the same factor. Finally, if 5 > 1 but the 6xi are of widely 
varying magnitudes, S is recalculated with only those results for which 5xi < 3N^^'^5x being 
retained; the recalculated scale factor is then applied in the usual way. In all three cases, no 
change is made to the original average x calculated with Eq. (jS)). 

The data for Qec include measurements of both individual Q^;^- values and the differences 
between pairs of Qbc- values. This required a two-step analysis procedure. We first treated 
the individual Qec- value measurements for each particular transition in the manner already 
described, obtaining an average result with uncertainty in each case, Xj ± 6xj, where the 
subscript j now designates a particular transition. For transitions unconnected by difference 
measurements, these uncertainties were scaled if necessary and then the values were quoted 
as final results. For those transitions involved in one or more difference measurements we 
combined their average Qec values, Xj±6xj, with the difference measurements, dk^Sdf^, in a 
single fitting procedure. If Mi is the number of transitions that are connected by difference 
measurements, and M2 is the number of those difference measurements, then we have a 
total of Ml + M2 input data values from which we need to extract a final set of Mi average 
Q _Bc- values, xj ± 6xj . We accomplish this by minimizing ? where 



with ji and j2 designating the two transitions whose Qec- value difference is determined in 
a particular d^ measurement. For each of these individual Q^c- values, we obtained its scale 
factor from Eq. where the used in that equation is now given by 




(6) 





and 
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where j is the particular transition of interest. The sum in i extends over all individual Qec- 
value measurements of transition j, and the sum in / extends over all doublet measurements 
that include transition j as one component. The resultant value of S was applied to the 
uncertainty, 6xj, with the same conventions as were described previously. 

B. Data tables 

The Q£;c- value data appear in Tables HI and HTl For the best-known nine superallowed 
decays - those of i°C, ^^O, ^SAP, ^^ci, ssj^m^ 42gp^ 46y^ so^n and ^^Co - the daughter 
nuclei are stable, and the most precise determinations of their Q^;^- values have come from 
direct measurements of that property via, for example, (p,n) or (^He,t) reactions. Such 
measurements are identified in column 3 of TableUby "Qscl-^fl)" and each individual result 
is itemized with its appropriate reference in the next three columns. The weighted average 
(see Eq. (jS}) of all measurements for a particular decay appears in column 7, with the cor- 
responding scale factor (see Eq. (jH))) in column 8. A few of these cases, like ^"^Cl and ^^V, 
have no further complications. There are other cases, however, in which Qec- value differ- 
ences have been measured in addition to the individual Q^jc- values. These measurements 
are presented in Table |n] They have been dealt with in combination with the direct Q ec- 
value measurements, as described in section HT^ (see, in particular, Eq. ((Zj)), with the final 
average Q^;^- value appearing in column 7 of Table H] and the average difference in column 4 
of Table |n] Both are flagged with footnotes to indicate the interconnection. 



TABLE I: Decay energies, QeC) for superallowed /3-decay branches. (See Table fVIIII for the correlation between the alphanu- 
meric reference code used in this table and the actual reference numbers.) 



Parent/Daughter Property^ 




Measured Energy, Qec (keV) 


Average value 




nuclei 


1 




2 3 


Energy (keV) 


scale 


= -1: 












lOC Qec{9s) 


3647.84 ± 0.34 


[Ba84] 


3647.95 ± 0.12 [Ba98] 


3647.94 ± 0.11 


1.0 




1740.15 ± 0.17 


[Aj88] 


1740.07 ± 0.02 2 


1740.07 ± 0.02 


1.0 


QEc(sa) 








1907.87 ± 0.11 




"O QEcigs) 


5143.35 ± 0.60 


[Bu61] 


5145.09 ± 0.46 [Ba62] 5145.57 ± 0.48 [Ro70] 








5142.71 ± 0.80 


[Vo77] 


5143.43 ± 0.37 [Wh77] 5144.34 ± 0.17 [To03] 


5144.29 ± 0.28 


2.1 
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TABLE I (continued) 



Parent/Daughter Property-'" Measured Energy (ke\') jYverase value 



nuclei 12 3 Energy (keV) scale 







fix(ciO+) 


2312.798 ± 0.011 [Aj91] 






2312.798 ± 0.011 








Qec (so) 










2831.18 ± 0.24 3 


2.5 


i^Ne 


18p 


ME(p) 


5316.8 ± 1.5 


[Ma94] 


5317.63 ± 0.36 [BI04b] 




5317.58 ± 0.35 


1.0 






ME{d) 


873.31 ± 0.94 


[Bo64] 


875.5 ± 2.2 [Ho64] 


876.5 ± 2.8 [Pr67] 












877.2 ± 3.0 


[Se73] 


873.96 ± 0.61 [Ro75] 




874.02 ± 0.48 


1.0 






Qec{9s) 


4438 ± 9 


[Pr63] 






4443.54 ± 0.60 


1.0 






E^{dO+) 


1 c\A 1 c^c; _|_ n nc 
iU4l:1.00 it U.Uo 


[iiyoj 






1041.55 ± 0.08 








Qec (so) 










3401.99 ± 0.60 




^^Mg 


22Na 


ME{p) 


-401.3 ± 3.0 


[Ha74c] 


-400.5 ± 1.3 4 




-400.6 ± 1.2 


1.0 






ME{d) 


-5184.3 ± 1.5 


[We68] 


-5182.5 ± 0.5 [Be68] 


-5181.3 ± 1.7 [An70] 












-5183.2 ± 1.0 


[Gi72] 


-5181.56 ± 0.16 [Mu04] 


-5181.08 ± 0.30 [Sa04] 


-5181.58 ± 0.27 


2.4 






Qec{9s) 


4781.64 ± 0.28 


[Mu04] 


4781.40 ± 0.67 [Sa04] 




4781.58 ± 0.25 


1.0 








657.00 ± 0.14 


[En98] 






657.00 ± 0.14 








Qec (so) 










4124.58 ± 0.29 




26 Si 


26 Al 


ME{p) 


-7159 ± 18 


[Mi67] 


-7149 ± 30 [Mc67] 


-7139 ± 30 [Ha68] 












-7145.5 ± 3.0 


[Ha74c] 






-7145.8 ± 2.9 


1.0 






ME{dO+) 


-11982.08 ± 0.26 


5 






-11982.08 ± 0.26 








Qec (so) 


4850 ± 13 


[Fr63] 






4836.9 ± 3.0 


1.0 


30g 


30p 


ME{p) 


-14060 ± 15 


[Mi67] 


-14054 ± 25 [Mc67] 


-14068 ± 30 [Ha68] 












-14063.4 ± 3.0 


[Ha74c] 






-14063.2 ± 2.9 


1.0 






ME{d) 


-20203 ± 3 


[Ha67] 


-20200.58 ± 0.40 [Re85] 




-20200.62 ± 0.40 


1.0 
















fii Q7 4 -1- 2 Q 










677.29 ± 0.07 [En98] 














QEcisa) 


5437 ± 17 


[Fr63] 






5459.5 =b 3.9 


1.3 


34 Ar 


34 Q 


IVI Hi \jj J 


-18380.2 ± 3.0 


[Ha74c] 


1 R'\77 1 n -1- n 41 riTpn9i 

-LOO ( / .-LU ^ U.^-L IJUcU^J 




1 ^7 -\- Ci do 


1.0 








-24440.01 ± 0.23 


6 






-94440 01 4- n 2S 


















6062 83 -1- 46 




38 Ca 


38k 


ME{p) 


-22056.0 ± 5.0 


[Se74] 






-22056.0 ± 5.0 








ME{dO+) 


-28670.20 ± 0.32 


5 






-28670.20 ± 0.32 








Qec (so) 










6614.2 ± 5.0 






42 Sc 


ME{p) 


-25121 ± 6 


[Mi67] 


-25086 ± 30 [Ha68] 


-25124 ± 13 [Zi72] 


-25120.7 ± 5.3 


1.0 






ME{d) 


-32121.55 ± 0.80 


5 






-32121.55 ± 0.80 








Qec (so) 










7000.9 ± 5.4 




T, = 0: 




















26 Mg 


Qec(9>^) 


4004.79 ± 0.55 


[De69] 


4004.41 ± 0.10 6 




4004.42 ± 0.10 


1.0 








228.:-!05 ± O.Oi:-! 


|Eu98] 






228.:-!05 ± O.Oi:-! 
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TABLE I (continued) 



Parent /D 


aughter Property^ 






Measured Energy (keV) 




Average value 




nuclei 




1 




2 


3 


Energy (keV) 


scale 






Qec {so) 


4232.71 ± 0.60 


fVo771 

[ V U 1 (J 


4232.19 ± 0.12 [Br94] 




4232.55 ± 0.17 ^ 


2.7 


34 


34 g 


Qec (sa) 


5490.3 ± 1.9 
5492.2 ± 0.4 


fRv73al 
[Vo77] 


5491.6 ± 2.3 [Ha74d] 
5491.65 ± 0.26 ^ 


5491.71 ± 0.54 [Ba77c] 


5491.78 ± 0.20 


1.0 




38 Ar 


QEcigs) 

E^{pO+) 
Qec (so) 


5914.76 ± 0.60 
130.4 ± 0.3 
6044.6 ± 1.5 


[Ja78] 
[En98] 
[Bu79] 


rrx A A oo 1 1 r> FTT r\ol 

6044.38 ± 0.12 [lla98J 




5914.76 ± 0.60 
130.4 ± 0.3 
6044.40 ± 0.11 


1.0 


42 Sc 


42 Ca 


Qec (so) 


6423.71 ± 0.40 


[Vo77] 


6425.84 ± 0.17 ° 




6425.63 ± 0.38 ' 


3.2 


46 V 


46 Ti 


Qec (sa) 


"7/^ COO 1 1 O 

7053.3 ± 1.8 


[Sq76] 


7050.41 ± 0.60 [Vo77] 




7050.71 ih 0.89 


1.6 


SOMn 


50 Cr 


Qec (sa) 


7632.8 ± 2.8 


[Ha74d] 


7631.91 ± 0.40 [Vo77] 




7632.43 ± 0.23 ^ 


1.0 


54Co 


54pg 


QEc(sa) 


8241.2 ± 1.8 


[Ho74] 


8245.6 ± 3.0 [Ha74d] 


8241.61 ± 0.60 [Vo77] 


8242.60 ± 0.29 ^ 


1.5 


62 Ga 


62 Zn 


Qec {so.) 


9171 ± 26 


[Da79] 






9171 ± 26 




66 As 


66 Ge 


Qec (so.) 


9550 ± 50 


[Da80] 






9550 ± 50 




™Br 


™Se 


Qec (so) 


9970 ± 170 


[Da80] 






9970 ± 170 




74Rb 


74Kr 


ME{p) 
ME{d) 
Qec (so.) 


-51905 ± 18 
-62330.3 ± 2.4 


[He02] 
[He02] 


-51915.2 ± 4.0 [Ke04] 
-62332.0 ± 2.1 [Ke04] 




-51914.7 ± 3.9 
-62331.3 ± 1.6 
10416.5 ± 4.2 


1.0 
1.0 



^ Abbreviations used in this column are as follows: "gs", transition between ground states; "sa", super- 
allowed transition; "p", parent; "d", daughter; "ME", mass excess; "Ex{0~^y\ excitation energy of the 
0+ (analog) state. Thus, for example, "Qec{scl)" signifies the Q^c-value for the superallowed transition, 
"ME{dy\ the mass excess of the daughter nucleus; and "ME{dO^), the mass excess of the daughter's 0+ 
state. 

^ Result based on references [Ba88] and [Ba89] . 

Average result includes the results of Qec pairs; see Table Hll 

Result based on references [Bi03] and [Se04] . 
^ Result obtained from data elsewhere in this table. 

6 Resuh based on references [IsSO], [A182], [Hu82], [Be85], [Pr90], [Ki91] and [Wa92]. 

Result based on references [Wa83] , [Ra83] and [Li94] . 
^ Result based on references [Zi87] and [Ki89]. 
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TABLE II: Qec-^^^^^ differences for superallowed /3-decay branches. Tliese data are also used 



as input to determine some of the average (5£;c-values Hsted in Tabled (See Table IVIIII for the 
correlation between the alphabetical reference code used in this table and the actual reference 
numbers.) 



T> J. 

Farent 


Parent 


(^EC2 — QeCI 


(keV) 


nucleus 1 


nucleus 2 


measurement 


average" 


140 

^^Al™ 

42SC 
42SC 

50Mn 


42SC 

5°Mn 

s^Co 

s^Co 


1401.68 ± 0.13 [Ko87] 

2193.5 ± 0.2 [Ko87] 

1207.6 ± 2.3 [Ha74d] 
1817.2 ± 0.2 [Ko87] 
610.09 ± 0.17 


1401.37 ± 0.29 
2193.09 ± 0.42 
1206.79 ± 0.44 
1816.97 ± 0.48 
610.18 ± 0.37 



"Average values include the results of direct Q^c-value measurements: see Tabled 



There are two cases, "^^p^^ and ^sj^*", in which the superallowed decay originates from an 
isomeric state. For both, there are Qsc- value measurements that correspond to the ground 
state as well as to the isomer. Obviously, these two sets of measurements are simply related 
to one another by the excitation energy of the isomeric state in the parent. In Table HI the 
set of measurements for the ground-state Q^c- value and for the excitation energy of the 
isomeric state appear in separate rows, each with its identifing property given in column 3 
and its weighted average appearing in column 7. In the row below, the average value given 
in column 7 for the superallowed transition is the weighted average not only of the direct 
superallowed Q^;^- value measurements in that row, but also of the result derived from the 
two preceeding rows. Note that in all cases the Qec- value for the superallowed transition 
appears in bold-face type. 

For those 11 superallowed decays that lead to radioactive daughter nuclei, there are very 
few direct measurements of the Qbc- value for the superallowed transition. In general, that 
Q_EC- value must be deduced from the measured mass excesses of the parent and daughter 
nuclei, together with the excitation energy of the analog 0^ state in the daughter. Each of 
these properties is identified in column 3 of Table U with the individual measurements of 
that property, their weighted average and a scale factor appearing in columns to the right. 



10 



TABLE III: Half-lives, ti/2) of superallowed /3-emitters. (See Table lVIIll for the correlation between 



the alphabetical reference code used in this table and the actual reference numbers.) 



Parent 
nucleus 






Measured half-lives, 4^/2 (ms) 




Average value 


1 




2 


3 


4 




scale 


T, = -1: 
















lOc 


19280 ± 20 


;Az74] 


19295 ± 15 [Ba90] 






19290 ± 12 


1.0 


14Q 


70480 ± 150 


;a172] 


70588 ± 28 [C173] 


70430 ± 180 [Az74] 


70684 ± 77 [Be78] 








70613 ± 25 


;Wi78] 


70560 ± 49 [GaOl] 


70641 ± 20 [Ba04] 




70616 ± 14 


1.1 


iSNe 


1690 ± 40 


;as7o] 


1670 ± 20 [A170] 


1669 ± 4 [A175] 


1687 ± 9 [Ha75] 


1672.1 ± 4.6 


1.3 


22Mg 


3857 ± 9 


;Ha75] 


3875.5 ± 1.2 [Ha03] 






3875.2 ± 2.4 


2.0 


26Si 


2210 ± 21 


;Ha75] 


2240 ± 10 [Wi80J 






2234 ± 12 


1.3 


30s 


1180 ± 40 


;Ba67] 


1220 ± 30 [Mo71] 


1178.3 ± 4.8 [Wi80] 




1179.4 ± 4.7 


1.0 


34 Ar 


844.5 ± 3.4 


;Ha74a] 


847.0 ± 3.7 [la03J 






845.6 ± 2.5 


1.0 


38Ca 


470 ± 20 


;Ka68] 


439 ± 12 [Ga69] 


450 ± 70 [Zju2,\ 


430 ± 12 [Wi80] 


440.0 ± 7.8 


1.2 


42Ti 


200 ± 20 


;Ni69] 


202 ± 5 [Ga69] 


1 11/1 FA ^^'f\'\ 

173 ± 14 [A169] 




198.8 ± 6.3 


1.4 


= 0: 

2G ^ym 


6346 ± 5 


[Fr69a] 


6346 ± 5 [Az75] 


6339.5 ± 4.5 [AI77] 


6346.2 ± 2.6 [Ko83] 


6345.0 ± 1.9 


1.0 




6345 ± 14 


[Sc05] 












34 Cl 


1526 ± 2 


[Ry73a] 


1525.2 ± 1.1 [Wi76] 


1527.7 ± 2.2 [Ko83] 


1527.1 ± 0.5 [Ia03] 


1526.77 ± 0.44 


1.0 


38 j^m 


925.6 ± 0.7 


[Sq75] 


922.3 ± 1.1 [Wi76] 


921.71 ± 0.65 [Wi78] 


924.15 ± 0.31 [Ko83] 








924.4 ± 0.6 


[BaOO] 


QOA A(\ -1- n 1 zl [Rtinf^l 
yZ4t.'4:D ^ U-l^ [lJ<iUOJ 
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Z.o 


42 Sc 


680.98 ± 0.62 


[Wi76] 


680.67 ± 0.28 [Ko97a] 






680.72 ± 0.26 


1.0 


46v 


422.47 ± 0.39 


[A177] 


422.28 ± 0.23 [Ba77a] 


422.57 ± 0.13 [Ko97a] 




422.50 ± 0.11 


1.0 


50Mn 


284.0 ± 0.4 


[Ha74b] 


282.8 ± 0.3 [Fr75] 


282.72 ± 0.26 [Wi76] 


283.29 ± 0.08 [Ko97a] 


283.24 ± 0.13 


1.8 


54 Co 


193.4 ± 0.4 


[Ha74b] 


193.0 ± 0.3 [Ho74] 


193.28 ± 0.18 [A177] 


193.28 ± 0.07 [Ko97a] 


193.271 ± 0.063 


1.0 


62 Ga 


115.95 ± 0.30 


[A178] 


116.34 ± 0.35 [Da79] 


115.84 ± 0.25 [Hy03] 


116.19 ± 0.04 [B104a] 


116.175 ± 0.038 


1.0 




116.09 ± 0.17 


[Ca05] 












66 As 


95.78 ± 0.39 


[A178] 


95.77 ± 0.28 [Bu88] 






95.77 ± 0.23 


1.0 


70Bi, 


80.2 ± 0.8 


[A178] 


78.54 ± 0.59 [Bu88] 






79.12 ± 0.79 


1.7 


74Rb 


64.90 ± 0.09 


[OiOl] 


64.761 ± 0.031 [BaOl] 






64.776 ± 0.043 


1.5 



The average Q£;c-value listed for the corresponding superallowed transition is obtained from 
these separate averages. If a direct measurement of the superallowed Qec- value exists, then 
it is also included in the final average. 

Especially in these latter 11 cases, it might be imagined that it would have been sufficient 



for us to use the 2003 Mass tables 



2Q| to derive the Qec- values of interest. There are, how- 



ever, significant differences in our approach. We have included all pertinent measurements 



for each property as described in section III AI typically, only a subset of the available data 
is included as input to the mass tables. Furthermore, we have examined each reference in 
detail and either accepted the result, updated it to modern calibration standards or rejected 
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TABLE IV: Branching ratios, R, for superaHowed /3-transitions. ( See Table IVlTll for the correlation 



between the alphabetical reference code used in this table and the actual reference numbers.) 



Parent /D; 


aughter 


Daughter state 


Measured Branching Ratio, R (%) 


Average value 


nuclei 


(MeV) 


1 


2 


R (%) scale 


T, = -1: 














lOg 


2.16 
1.74 


0+0.0008 [(.^72] 

1.468 ± 0.014 [Ro72] 
1.465 ± 0.009 [Kr91] 
1.4665 ± 0.0038 [Fu99] 


1.473 ± 0.007 [Na91] 
1.4625 ± 0.0025 [Sa95] 


r, + 0.0008 

—0 

1.4646 ± 0.0019 1.0 


14q 


14N 


gs 

3.95 
2.31 


0.60 ± 0.10 [Sh55] 
0.61 ± 0.01 [Si66] 
0.062 ± 0.007 [Ka69] 
0.053 ± 0.002 [Hc81] 


0.65 ± 0.05 [Fr63] 
0.058 ± 0.004 [Wi80] 


0.611 ± 0.010 1.0 

0.0545 ± 0.0019 1.1 
99.334 ± 0.010 


i8Ne 


18 p 


1.04 


9 ± 3 [Fr63] 


7.70 ± 0.21" [Ha75] 


7.70 ± 0.21 1.0 


22Mg 


22Na 


0.66 


54.0 ± 1.1 [Ha75] 


53.15 ± 0.12 [Ha03] 


53.16 ± 0.12 1.0 


26 Si 


26A1 


1.06 
0.23 


21.8 ± 0.8 [Ha75] 




21.8 ± 0.8 

75.09 ± 0.92" 


30 c 


30 p 


gs 

0.68 


20 ± 1 [Fr63] 




20 ± 1 
77.4 ± 1.0" 


34 A „ 


34pi 


0.67 
gs 


o ACi _L n 1 n '7/1 f^l 
z.4y ih U.IU [rla/4aj 




2.49 ± 0.10 
94.45 ± 0.25" 


42Ti 


42 So 


0.61 
gs 


J- 1 A FA IfiOl 
OD it 14 [/vioyj 




56 ± 14 
43 ± 14" 


= 0: 














26 Mg 


gs 


>99.997 [Ki91] 




100.000lOoo3 


34 Cl 


34g 


gs 


>99.988 [Dr75] 




100.000l°oi2 


38 p^m 


38 Ar 


3.38 
gs 


<0.0019 [Ha94] 
>99.998 




n+0.002 

100.000t°.oo2 


42 Sc 


42 Ca 


1.84 
gs 


0.0063 ± 0.0026 [In77] 
0.0103 ± 0.0031 [Sa80] 


0.0022 ± 0.0017 [De78] 
0.0070 ± 0.0012 [Da85] 


0.0059 ± 0.0014 1.6 
99.9941 ± 0.0014 


46y 


46 Ti 


2.61 
4.32 

SGT'' 
gs 


0.0039 ± 0.0004 [Ha94] 
0.0113 ± 0.0012 [Ha94] 
<0.004 




0.0039 ± 0.0004 
0.0113 ± 0.0012 

r, + 0.004 

"-0 

99 9848+" '"'^^ 


SOMn 


50cr 


3.63 
3.85 
5.00 
gs 


0.057 ± 0.003 [Ha94] 
<0.0003 [Ha94] 
0.0007 ± 0.0001 [Ha94] 




0.057 ± 0.003 

r, + 0.0003 

0.0007 ± 0.0001 
99.9423 ± 0.0030 


54Co 


54pg 


2.56 

SGT*" 
gs 


0.0045 ± 0.0006 [Ha94] 
<0.03 




0.0045 ± 0.0006 

r, + 0.03 

"-0 

99 9955+" '"""' 


62 Ga 


62Zn 


SGT*" 
gs 


0.15l»;J^ [Hy03],[B102] 






-i-n n.^i 
99.85^111 


74Rb 


74Kr 


SGT'' 
gs 


0.50 ± 0.10 [Pi03] 




0.50 ± 0.10 
99.50 ± 0.10 



"^Result also incorporates data from Table IVl 

''designates total Gamow- Teller transitions to levels not explicitly listed; values were derived with the help 



of calculations in [Ha02l. 
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TABLE V: Relative intensities of /3-delayed 7-rays in the superallowed /?-decay daughters. These 
data are used to determine some of the branching ratios presented in Table HVl (See Table rvTTTI for 
the correlation between the alphabetical reference code used in this table and the actual reference 



numbers.) 



Parent /Daughter 
nuclei 


daughter 

ratios" 


Measured 


7-ray Ratio 


Average value 




1 


2 


Ratio 


scale 


l8Ne 


18p 


7660/71042 


0.021 ± 0.003 [Ha75] 


0.0169 ± 0.0004 [Hc82] 












0.0172 ± 0.0005 [Ad83] 




0.0171 ± 0.0003 


1.0 




26 Al 


71622 /7829 


0.149 ± 0.016 [Mo71] 


0.134 ± 0.005 [Ha75] 












0.1245 ± 0.0023 [Wi80] 




0.1265 ± 0.0036 


1.7 






71655 /7829 


0.00145 ± 0.00032 [WiSO] 




0.0015 ± 0.0003 








71843/7829 


0.013 ± 0.003 [Mo71] 


0.016 ± 0.003 [Ha75] 












0.01179 ± 0.00027 [WiSO] 




O.OllS ± 0.0003 


1.0 






"7951 9 / 'vs^q 


0.00282 ± 0.00010 [Wi80] 




0.002s ± 0.0001 








7total/7829 






0.1426 ± 0.0036 




30 s 


30p 


7709/7677 


0.006 ± 0.003 [Mo71] 


0.0037 ± 0.0009 [WiSO] 


0.0039 ± 0.0009 


1.0 






72341 /7677 


0.033 ± 0.002 [Mo71] 


0.0290 ± 0.0006 [WiSO] 


0.0293 ± 0.0011 


1.9 






73019/7677 


0.00013 ± 0.00006 [Wi80] 




0.0001 ± 0.0001 








7total/7677 






0.0334 ± 0.0014 




34 


34g 


7461 /7666 


0.28 ± 0.16 [Mo71] 


0.365 ± 0.036 [Ha74a] 


0.361 ± 0.035 


1.0 






72580 /7666 


0.38 ± 0.09 [Mo71] 


0.345 ± 0.01 [Ha74a] 


0.345 ± 0.010 


1.0 






73129/7666 


0.67 ± 0.08 [Mo71] 


0.521 ± 0.012 [Ha74a] 


0.524 ± 0.022 


l.S 






7total/7666 






1.231 ± 0.043 




42Ti 


42 Sc 


72223 /7611 


0.012 ± 0.004 [Ga69] 




0.012 ± 0.004 








7total/7611 


2 X 0.012 [En90] 




0.024 ± 0.008 





''7-ray intensities are denoted by 7^;, where E is the 7-ray energy in keV. 

it for cause. The updating procedures are outlined, reference by reference, in Table EH and 
the rejected results are similarly documented in Table rVTTl With a comparatively small data 
set, we could afford to pay the kind of individual attention that is impossible when one is 
considering all nuclear masses. 

The half-life data appear in Table IIIII in similar format to Table lU For obvious reasons, 
half-life measurements do not lend themselves to being updated. Consequently, a number 
of mostly pre-1970 measurements have been rejected because they were not analyzed with 
the "maximum-likelihood" method. The importance of using this technique for precision 
measurements was not recognized until that time and, without access to the primary 
data, there is no way a new analysis can be applied retroactively. All rejected half-life 
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TABLE VI: References for which the original decay-energy results have been updated to incorporate 
the most recent calibration standards. (See Table Ivnil for the correlation between the alphabetical 
reference code used in this table and the actual reference numbers.) 

References (parent nucleus)" Update procedure 

• We have converted all original (p,n) threshold measurements to Q-values 
using the most recent mass excesses [Au03]. 

• These (p,n) threshold measurements have been adjusted to reflect recent 
calibration a-energies [Ry91] before being converted to Q-values. 

• Before conversion to a Q-value, this (p,n) threshold was adjusted to reflect a 
new value for the '^Li(p,n) threshold [Wh85], which was used as calibration. 

• This (p,n) threshold was measured relative to those for ^'^C and ^*0; we have 
adjusted it based on average Q-values obtained for those decays in this work. 

• Before conversion to a Q-valuc, this (p,n) threshold was adjusted to reflect the 
modern value for the ^^Cl(p,n) threshold [Au03], which was used as calibration. 

• These ^^C('^He,n) threshold measurements have been adjusted for updated 
calibration reactions based on current mass excesses [Au03]. 

• These ('^He,t) reaction Q-values were calibrated by the ^'^Al(''He,t) reaction 
to excited states in ■^'^Si; they have been revised according to modern mass 
excesses [AuOS] and excited-state energies [En98]. 

• This *-'^Ca(p,7) reaction Q-valuc was measured relative to that for *''Ca(p,7); 
we have slightly revised the result based on modern mass excesses [Au03]. 

• These (p,t) reaction Q-values have been adjusted to reflect the current Q- 
value for the ^^0{p,t) reaction [AuOS], against which they were calibrated. 

"These references all appear in Tabled under the appropriate parent nucleus. 

measurements are also documented in Table IVlTl 

Finally, the branching-ratio measurements are presented in Table IIVI The decays of the 
Tj = parents are the most straightforward since, in all these cases, the superallowed branch 
accounts for >99.5% of the total decay strength. Thus, even imprecise measurements of the 
weak non-superallowed branches can be subtracted from 100% to yield the superallowed 
branching ratio with good precision. For the higher- Z parents of this type, part icularly ^^Ga 
and heavier, it has been shown theoretically j8l| and experimentally [90|, Ill7l | that numer- 
ous very-weak Gamow- Teller transitions occur, which, in total, can carry significant decay 
strength. Where such unobserved transitions are expected to exist, we have used a combi- 
nation of experiment and theory to account for the unobserved strength, with uncertainties 
being adjusted accordingly. 



• Bo64(i8Ne), Ba84(ifC), Br94(26Al'") 
Ba98 C°C), Ha98 (^SK""), To03 C'^O) 

• Ry73a(34ci), Ho74(5*Co), Sq76 (^SV) 
Ba77c(34ci), Wh77(i4o) 

• Pr67(i8Ne) 

• Ja78 (38 K") 

• Bu79 (38k™) 

• Bu61 ("O), Ba62("0) 

• Ha74d (34C1, ^OMn, ^^Co) 

• Ki89(*2Sc) 

• Ha74c (22Mg, 26Si, 30g^ 34Ar)^ 
Se74(38Ca) 
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TABLE VII: References from which some or ah results have been rejected. (See Table IVIIII for 
the correlation between the alphabetical reference code used in this table and the actual reference 
numbers.) 

References (parent nucleus) Reason for rejection 

• No calibration is given for the measured (p,t) reaction Q-values; update 
is clearly required but none is possible. 

• Calibration reaction Q-values have changed but calibration process is too 
complex to update. 

• P.H. Barker (co-author) later considered that inadequate attention had 
been paid to target surface purity [Ba84]. 

• P.H. Barker (co-author) later stated [Ba84] that the (p,t) reaction Q-value 
could not be updated to incorporate modern calibration standards. 

• The result in [WhSl] was updated in [Ba98] but then eventually virithdrawn 
by P.H. Barker (co-author) in [To03]. 

• Quoted uncertainties are too small, and results likely biased, in light of 
statistical difficulties more recently understood (see [Fr69a]). In particular, 
"maximum-likelihood" analysis was not used. 

• All four quoted half-lives are systematically higher than more recent and 
accurate measurements. 

• P.H. Barker (co-author) later considered that pile-up had been 
inadequately accounted for [Ba90]. 

• "Maximum-likelihood" analysis was not used. 

• Numerous impurities present; result is obviously wrong. 

The branching ratios for decays from = — 1 parents are much more challenging to 
determine, since the superallowed branch is usually one of several strong branches - with 
the notable exception of ^^O - and, in two of the measured cases, it actually has a branching 
ratio of less than 10%. The decays of ^'^F, ^^Si, ^°S, ^^Ar and ^^Ti required special treatment 
in our presentation. In each of these five cases, the absolute branching ratio for a single 
/5-transition has been measured. The branching ratios for other /^-transitions must then be 
determined from the relative intensities of /3-delayed 7-rays in the daughter. The relevant 
7-ray intensity measurements appear in Table 1^ with their averages then being used to 
determine the superallowed branching-ratio averages shown in bold type in Table IIVI These 



1. Decay-energies: 

• Pa72(30s, 38 Ca) 

• No74 (22Mg) 

• Ro74(iOC) 

• Ba77b(i''C) 

• WhSl and Ba98 C^O) 

2. Half-lives: 

• Ja60(26Al'"), He61(i*0), Ba62(i4o)^ 
Ea62(iOC), Ba63(i0c), Fr63(i4o, 2631), 
FreSbC'^Sc, 4«V, 50Mn), Si72("0) 

• Ha72a(26Al™, 34ci, SS-^m^ 42g^) 

• Ro74(iOC) 

• Ch84(38K'") 

3. Branching-ratios: 

• Fr63(26s) 
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TABLE VIII: Reference key, relating alphabetical reference codes used in Tables Hillll to the actual 
reference numbers. 



Table 


Reference 


Table 


Reference 


Table 


Reference 


Table 


Reference 


Table 


Reference 


Table 


Reference 


code 


number 


code 


number 


code 


number 


code 


number 


code 


number 


code 


number 


Ad83 




Ba98 


[33] 


Dr75 


[58] 


Hc61 


[83] 


Mc67 


[1081 


Se74 


[133] 


Aj88 


m 


BaOO 


[34] 


Ea62 


[59] 


He81 


[84] 


Mi67 


[109] 


Se04 


[1341 


Aj91 


[10] 


Bad 


[35] 


En90 


[60] 


He82 


[85] 


Mo71 


[110] 


Sh55 


[1351 


A169 




Ba04 


[36] 


En98 


[61J 


He02 


[86] 


Mu04 


[111] 


Si66 


[1361 


A170 


[12] 


Ba05 


[^ 


Fr63 


m 


Ho64 


m 


Na91 


[112] 


Si72 


[1371 


A172 


yj 


Be68 


[22] 


Fr65b 


m 


Ho74 


m 


Ni69 


[112] 


Sq75 


[1381 


A175 


[14] 


Be78 


[^ 


Fr69a 


m 


Hu82 


m 


No74 


[114] 


Sq76 


[1391 


A177 


[15] 


Be85 


[40] 


Pr75 


[65] 


Hy03 


[90] 


OiOl 


[115] 


Ti95 


[1401 


A178 


[16] 


Bi03 


[41] 


Pu99 


[66] 


Ia03 


[91] 


Pa72 


[116] 


To03 


[1411 


A182 


[17] 


B102 


[42] 


Ga69 


[67] 


In77 


[92] 


P103 


[117] 


Vo77 


[1421 


An70 


[IJ 


B104a 


[42] 


GaOl 


[£2] 


Is80 


m 


Pr67 


[112] 


Wa83 


[1431 


As70 


[19] 


B104b 


[44] 


Go72 


[70] 


Ja60 


[94] 


Pr90 


[119] 


Wa92 


[1441 


Au03 


[20] 


Bo64 


[45] 


Gi72 


[69] 


Ja78 


[95] 


Ra83 


[120] 


We68 


[1451 


Az74 


[21] 


Br94 


[46] 


Ha67 


[71] 


Ka68 


[96] 


Re85 


[121J 


Wh77 


[1461 


Az75 


[22] 


Bu61 


[47] 


Ha68 


[72] 


Ka69 


[97] 


Ro70 


[122] 


Wh81 


[1471 


Ba62 


[23] 


Bu79 


[48] 


Ha72a 


[73] 


Ke04 


[98] 


Ro72 


[123] 


Wh85 


[1481 


Ba63 


[24] 


Bu88 


[49] 


Ha74a 


[74] 


Ki89 


[99] 


Ro74 


[124] 


Wi76 


[1491 


Ba67 


[25] 


Ca05 


[50] 


Ha74b 


[75] 


Ki91 


[100] 


Ro75 


[125] 


Wi78 


[1501 


Ba77a 


[26] 


Ch84 


[51] 


Ha74c 


[76] 


Ko83 


[101] 


Ry73a 


[126] 


Wi80 


[1511 


Ba77b 


[27] 


C173 


[52] 


Ha74d 


[77] 


Ko87 


[102] 


Ry91 


[127] 


Zi72 


[1521 


Ba77c 


[28] 


Da79 


[53] 


Ha75 


[78] 


Ko97a 


[103] 


SaSO 


[128] 


Zi87 


[1531 


Ba84 


[29] 


Da80 


[54] 


Ha94 


[79] 


Ko97b 


[104] 


Sa95 


[129] 






Ba88 


[30] 


Da85 


[55] 


Ha98 


[80] 


Kr91 


[105] 


Sa04 


[130] 






Ba89 


[31] 


De69 


[56] 


Ha02 


[81] 


Li94 


[106] 


Sc05 


[131] 






Ba90 


[32] 


De78 


[57] 


Ha03 


[82] 


Ma94 


[107] 


Se73 









cases are also flagged with a table footnote. 

III. THE Tt VALUES 

Having surveyed the primary experimental data, we now turn to producing a set of ft- 
values for the twenty superallowed transitions being considered. The statistical rate function, 
/, for each transition depends primarily on the charge of the daughter nucleus, Z, and on 
the Q^;^- value to the fifth power. Consequently the uncertainty in the value of / due to 
the experimental uncertainty in Qec is given by (A///) ^ 5{AQec/Qec)- Our goal in 
computing / therefore is to ensure that the computation itself yields percentage errors much 
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TABLE IX: Derived results for superallowed Fermi beta decays. 



Parent 




Pec 


Partial half-life 










nucleus 


/ 


(%) 


t{ms) 




S'r (%) 


Sc - Sns (%) 




= -1: 


















2.3009 ±0.0012 


0.297 


1321000 ± 1900 


3039.5 ±4.7 


1.652 ±0.004 


0.540 ± 0.039 


3073.0 ±4.9 


14 o 


42.772 ± 0.024 


0.088 


71151 ± 16 


3043.3 ± 1.9 


1.520 ±0.008 


0.570 ± 0.056 


3071.9 ±2.6 


iSNe 


134.48 ±0.15 


0.081 


21730 ± 590 


2922 ± 80 


1.484 ±0.012 


0.910 ±0.047 


2938 ± 80 


22 Mg 


418.44 ±0.18 


0.069 


7295 ± 17 


3052.4 ± 7.2 


1.446 ±0.017 


0.505 ±0.024 


3080.9 ±7.4 


26Si 


1023.3 ±3.7 


0.064 


2978 ± 40 


3047 ± 42 


1.420 ±0.023 


0.600 ± 0.024 


3072 ± 42 


30s 


1967.1 ±7.8 


0.066 


1524 ± 21 


2998 ± 44 


1.405 ± 0.029 


1.125 ±0.039 


3006 ± 44 


34 Ar 


3414.2 ± 1.5 


0.069 


896.0 ±3.5 


3059 ± 12 


1.394 ±0.035 


0.825 ± 0.044 


3076 ± 12 


38Ca 


5338 ± 22 


0.075 






1.397 ±0.042 


0.910 ±0.053 




42Ti 


7043 ± 30 


0.088 


470 ± 160 


3300 ±1100 


1.412 ±0.050 


1.015 ±0.110 


3300 ± 1100 


= 0: 
















26 j\^jm 


478.20 ±0.11 


0.082 


6350.2 ± 1.9 


3036.7 ± 1.2 


1.458 ±0.020 


0.261 ±0.024 


3072.9 ± 1.5 


34 CI 


1996.39 ±0.41 


0.080 


1527.99lo;l7 


3050.5 ± 1.1 


1.425 ±0.032 


0.720 ±0.039 


3071.7 ±1.9 




3298.10 ±0.33 


0.085 


925.11 ±0.27 


3051.1 ± 1.0 


1.423 ±0.039 


0.720 ±0.047 


3072.2 ±2.1 


42 Sc 


4470.03 ± 1.46 


0.099 


681.43 ±0.26 


3046.0 ± 1.5 


1.437 ± 0.047 


0.460 ± 0.047 


3075.6 ±2.5 


46v 


7200.0 ± 5.0 


0.101 


422.99 ±0.11 


3045.5 ± 2.2 


1.429 ±0.054 


0.465 ±0.033 


3074.7 ±3.0 


SOMn 


10731.2 ± 1.8 


0.107 


283.71 ±0.13 


3044.5 ±1.5 


1.429 ±0.062 


0.547 ±0.037 


3071.1 ±2.7 


e2Ga 
««As 


1 ^7AQ -t- n 
26250 ± 400 
31610 ±890 


111 
0.137 
0.156 


iqq /inK+0.063 

116.509in?^ 


3058 ± 47 


1.428 ±0.071 
1.445 ±0.087 
1.457 ± 0.095 


0.639 ± 0.043 
1.42 ±0.16 
1.45 ±0.16 


Qn71 9 -t- 9 
3058 ± 47 


70Br 


38600 ± 3600 


0.175 






1.47 ±0.11 


1.41 ±0.20 




74Rb 


47280 ± 100 


0.194 


65.227 ±0.078 


3083.8 ± 7.5 


1.49 ±0.12 


1.50 ±0.41 


3083 ± 15 












Average (best 12), J^t 


3072.7 ±0.8 
















0.42 



less than those due to the uncertainty in the Q^c- value, which can be <0.02% in the best 
cases. To this end we have written a new code, the details of which are given in Appendix 
1X1 Our final /-values and their uncertainties are recorded in column two of Table IIXI 

The partial half-life, t, for each transition is obtained from its total half-life, ti/2, and 
branching ratio, R, according to the formula 

t = ^ (1 + Pec) (9) 
where Pec is the calculated electron-capture fraction. The evaluation of Pec is discussed 



17 



by Bambynek et al. 154| and is based on the equation 



-EC = \n 



Y.l3l{WEc-\W^\fB, 



If (10) 



The sum extends over all atomic subshells from which an electron can be captured. The 
factor I3x is the Coulomb amplitude of the appropriate bound-state electron radial wave 
function; Wec is the Q^;^- value expressed in electron rest-mass units; Wx is the x-subshell 
binding energy also in electron rest-mass units; and B.^ takes account of the effects of electron 
exchange and overlap. We have computed Pec for the cases of interest here using our Qec- 
values from Table U and the values of fi'^Bx and Wx from, respectively, Tables 1 and 2 of 
Appendix F in ref. 155[. The Pec results are shown (as percentages) in column three of 
Table HXl Based on experimental tests of such Pec calculations |l54| . we expect these results 
to be accurate to a few parts in 100; thus they do not contribute perceptibly to the overall 
uncertainties. Partial half-lives derived from Eq. 0, and corresponding ft values appear in 
columns four and five. 

To obtain jFt-values according to Eq. Q we must now deal with the small correction 
terms. The term 5^ has been calculate d from st andard QED, and is currently evaluated to 



order Za^ and estimated in order Z'^a^ 



, 157j |: its values, listed in column six of Table HXt 
are around 1.4% and can be considered to be very reliable. The structure-dependent terms 
5ns and 5c, have also been calculated in the past but at various times over three decades 
and with a variety of different models. Th eir u ncertainties are larger. This topic has been 
reviewed recently by Towner and Hardy jlSS^ . who presented new calculations of these 
corrections in which consistent model spaces and approximations have been used for both 
correction terms. The results of these new calculations are recorded in column seven of 
Table HXl Finally, the resulting jFt- values are listed in column eight. 

A. CVC test 

We are now ready to test the CVC assertion that the Tt values should be constant for all 
nuclear superallowed transitions of this type. The data in Table HXl clearlv satisfy the test; 
the weighted average of the 12 most-precise results (with "statistical" uncertainty only) is 

= 3072.7 ± 0.8s (11) 
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FIG. 1: J-t values plotted as a function of the charge on the daughter nucleus, Z. The shaded 
horizontal band gives one standard deviation around the average J-t value, Eq. 

with a corresponding chi-square per degree of freedom of x^/^ = 0.42. In Fig. we plot the 
same 12 values, all of whose statistical accuracy is better than 0.5%. It is evident from both 
the figure and the table that the data form a consistent set, thus verifying the expectation 
of the CVC hypothesis at the level of 3 x 10~^, which is the fractional uncertainty quoted in 
Eq. (Illj). This is a 30% improvement over the results from our last survey in 1990 6] and 
is principally due to improvements in the experimental data themselves. 



B. .Ft- value error budgets 

In Figs. 121 and 01 we show the contributing factors to the individual JFt- value uncertainties. 
For the most precise data, ^^O to ^^Co, which appear in Fig. |2l the theoretical uncertainties 
are greater than, or comparable to, the experimental ones. The nuclear-structure-dependent 
correction, 6c — Sns, contributes an almost constant uncertainty of 4 parts in 10^ across 
these nuclei, while the nucleus-dependent radiative correction, 5^, has an uncertainty that 
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FIG. 2: Summary histogram of the fractional uncertainties attributable to each experimental and 
theoretical input factor that contributes to the final J^t values for the nine most precise superallowed 
decay data. 

grows as Z^. This is becaus e the contribution to 5p from order Z'^a^ has only been estimated 

n 

from its leading logarithm jl5fi| and the magnitude of this estimate has been taken as the 
uncertainty in 5^. In fact, for ^''Mn and ^^Co this becomes the leading uncertainty, indicating 
that a closer look at the order Z'^a'^ contribution would now be worthwhile. For the eight 
precise data, the experimental branching ratios are > 99% and have very small associated 
uncertainties with the exception of ^^Co, which has a 3 x 10^ fractional uncertainty attributed 
to its branching ratio. This is because ^^Co is predicted to have several weak Gamow- Teller 
branches that have not yet been observed. We have used an estimate of the strength of 

n 

the missing branches, taken from a shell-model calculation [81], to assign an uncertainty 
to the superallowed branching ratio. Missing weak branches become a larger issue for the 
heavier-mass nuclei with A > 62, where they contribute significantly to the branching-ratio 
uncertainty. 
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FIG. 3: Summary histogram of the fractional uncertainties attributable to each experimental and 
theoretical input factor that contributes to the final J^t values for the twelve other superallowed 
decay data. Where the error is shown as exceeding 60 parts in 10^, no useful experimental mea- 
surement has been made. 

For the less precisely known decay of ^''C, and for the twelve decays depicted in Fig. 01 



1 th e single exception of 
158| . resulting in a larger 



the predominant uncertainties are all experimental in origin wit 
^^Rb, for which the nuclear-structure calculation is quite difficult 
uncertainty on 6c — Sns- Many of the experimental Q-values, half-lives and branching 
ratios have yet to be measured for the cases in Fig. El but recent advances in experimental 
techniques are likely to change this situation dramatically within the next few years. 



21 



C. Accounting for systematic uncertainties 



So far, we have dealt with the inter-nuclear behavior of ^t-values, examining their con- 
stancy as a test of CVC With that test passed at high precision, we are now in a position to 
use the average JFt-value obtained from these concordant nuclear data to go beyond nuclei, 
obtaining first the vector coupling constant (see Eq. Q) and then the Vud matrix element. 
Before doing so, however, we must address one more possible source of uncertainty. Though 
the average J-'t value given in Eq. (Illj) includes a full assessment of the uncertainties at- 
tributable to experiment and to the particular calculations used to obtain the correction 
terms, it does not incorporate any provision for a common systematic error that could arise 
from the type of calculation chosen to model the nuclear-structure effects. In this section 
we look more critically at the nuclear-structure-dependent corrections, and in particular at 
the isospin-symmetry-breaking correction^, Sc- 

There have been a number of previous calculations of 6c besides those of ours 



Hartree-Fock cal cula tions of Ormand and Brown 
Giai and Suzuki 



1^ 



, R PA calculations of Sagawa, van 
16l|, and Woods-Saxon calculations 



160 1. R-matrix calculations of Barker 
i 

of Wilkinson [162], to name some of the more recent publications. Of these, we will only 
retain the Ormand-Brown (OB) calculations since they, like ours (TH), are constrained to 
reproduce other isospin properties of the nuclei involved: They reproduce the measured 
coefficients of the relevant isobaric multiplet mass equation, the known proton and neutron 
separation energies, and the measured ft values of weak non-analog 0^ 0^ transitions 
79|, where they are known. The other calculations are not constrained by experiment in 
any way and thus offer no independent means to assess their efficacy. 

Unfortunately, calculations of 6c by OB are not available for all the cases listed in Ta- 
ble |IXJ so we must concentrate on the nine most precise data: ^'^C, ^^O, ^^Al™, '^^Cl, ^^K^, 
^^Sc, ^^V, ^°Mn, and ^"^Co. When the TH values of 6c are used, the average J-'t-value for 
these nine cases alone is J-'t = 3072.6 ± 0.8 s with x^/^ = 0.35. When OB values are used 
for 6c instead, the weighted average is J-'t = 3074.5 ± 0.8 s with x^/z/ = 0.92. Although 
the chi-square with the OB values is nearly a factor of three worse, we do not argue that 



^ The reason we do not consider further the nuclear-structure-dependent radia tive correction, Sns, is that 
it is very small for the series of transitions that have Tz — parent states [ifisj . Of the nine precisely 
known transitions we are concentrating on, seven are of this type. 



22 



this is sufficient reason to reject the OB calculation. Rather, we observe that the OB values 
of 5c are systematically smaller and hence the jFt values systematically larger than ours. 
Evidently there is a systematic difference between our Woods-Saxon and OB's Hartree-Fock 
calculations of 5c and that difference should be accounted for in the final result. Thus, we 
adopt the average of these two results for our recommended jFt, and assign a systematic 
uncertainty equal to half the spread between them: viz 

Tt = 3073.5 ± O.Sstat ± 0.9syst s 

= 3073.5 ± 1.2 s, (12) 

where the two errors have been combined in quadrature. 



IV. IMPACT ON WEAK-INTERACTION PHYSICS 



A. The Value of 

With a mutually consistent set of J-'t values, we can now use their average value in 
Eq. (fT^ to determine the vector coupling constant, Gy, from Eq. The value of Gy itself 
is of little interest, but it can be related to the weak interaction constant for the purely 
leptonic muon decay, Gp, to yield the much more interesting up-down matrix element of the 
Cabibbo-Kobayashi-Maskawa (CKM) quark-mixing matrix^: Gv = G^Vud- The relation we 
use is 

K 

2GUl + Al)J^t' ^ ^ 

where is the nucleus-independent radiative correction. The currently accepted value for 



this correction is derived from the expression 163l Il64 1 



K = 7^[^ ln(mz/mp) + \ia{mp/m^) + 2CBorn] + ■ • • , (14) 

where the ellipses represent further small terms of order 0.1%. Here is the Z-boson 
mass, rup the proton mass, rriA the mass parameter in the dipole form of the axial-vector 



^ More completely we could write Gv — GpVudgv{q^ 0), where gv is the vector form factor given in 
Eq. (|A18]I . or as Gv = GpVudCv, where Cy is the vector coupling constant in the Jackson, Treiman and 
Wyld [163 Hamiltonian in Eq. with g^iq'^ ^ 0) = Gy = 1. 
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form factor, and CBorn is the universal order-a axial-vector contribution. The various terms 
have the values 

= 2.12 - 0.03 + 0.20 + 0.10%, (15) 

with the first term, the leadin g log arithm, being essentially unambiguous in value. The final 
value recommended by Sirlin (l64| is 

= (2.40 + 0.08)%. (16) 

The uncertainty is almost entirely due to the value selected for the axial-vector form factor 
mass, which Sirlin argues should lie in the range (ma^/2) < rrij^ < 2ma^, where is the 
physical oi meson mass. 

Using the Particle Data Group (PDG) [3] value for the weak interaction coupling constant 
from muon decay of G^/{hcf = (1.16639 + 0.00001) x 10"^ GeV-^, we obtain from Eq. (HHj) 
the result 

iK^r = 0.9482 + 0.0008. (17) 

Note that the total uncertainty here - 0.00083, if the next significant figure is included - is 
almost entirely due to the uncertainties contributed by the theoretical corrections. By far 
the largest contribution, 0.00074, arises from the uncertainty in A^; 0.00031 comes from the 
nuclear-structure-dependent corrections 6c — Sns (principally from the systematic difference 
between the OB and TH calculations discussed in Sect. IIIlT^ and 0.00012 is attributable to 
Only 0.00016 can be considered to be experimental in origin. 
The corresponding value of Vud is 

iKdl = 0.9738 + 0.0004, (18) 

a result that differs by two units in the last quoted digit from our previously recommended 



result 



This shift, well within the quoted one standard deviation, is due to the improve- 
ments in the experimental data and to our re-computing of the statistical rate function (see 
Appendix , in which a number of different parameter choices were made for the charge- 
density distribution, the oscillator length parameter for nuclear radial functions, and for the 
screening correction. Coincidentally, the value of Vud quoted in Eq. ()18j) is identical to the 
currently recommended PDG value [7], although our uncertainty is one digit smaller. 
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B. Unitarity of the CKM matrix 



The CKM matrix yields the transformation equations for a change of basis from quark 
weak-interaction eigenstates to quark mass eigenstates. As such, the CKM matrix must 
be unitary in order that the bases remain orthonormal. With the CKM matrix elements 
determined from experimental data, one important test they should satisfy is that they yield 
a unitary matrix. Currently, the sum of the squares of the top-row elements, which should 
equal one, constitutes the most demanding available test. With our experimental value for 
iKdP given in Eq. ^ and the PDG's recommended values [ll of |K.| = 0.2200 ± 0.0026 
and iKfel = 0.00367 ± 0.00047, this unitarity test yields: 

iK^r + iKsl' + iK^r = 0.9966 ± 0.0014 (19) 

The test fails by 2.4 standard deviations. Two-thirds of the assigned error is attributed to 
the uncertainty in |\4sp, viz. 0.0011, and one-third to the error in |\4dp, viz. 0.0008. The 
latter, as we have already demonstrated, is not predominantly experimental in origin, but 
is dominated by the uncertainty in the nucleus-independent radiative correction, A^. 

A recent measurement of the —>■ jr^e^Uf. {K^^) branching ratio from the Brookhaven 
E865 experiment P obtains Vus = 0.2272 ± 0.0030. Although this result is included in the 
PDG average value, it is considerably higher than the older experimental results from K^^ 
and decays, with which it is inconsistent. Experiments now in progress should help 
clarify the situation. If, for the moment, we adopt the E865 value for Vus rather than the 
PDG average, then the result in Eq. ()19|) is modified to 

iK^r + \Vus\' + \Vub\^ = 0.9999 ± 0.0016 (20) 

and unitarity is fully satisfied. 

Another, at present, less demanding test is to examine the first column of the CKM 
matrix. The PDG value for Vd is 0.224 ± 0.012, but little is known about Vtd other than it 
is expected to lie in the range 0.0048 < Vtd < 0.014. In this range it has negligible impact 
on the unitarity sum. With our value of |V^dP, this unitarity sum becomes 

\Vud\^ + \Vcd\^ + \Vtd\^ = 0.9985 ± 0.0054 (21) 

Here the error is given entirely by the uncertainty in the value of Vcd and unitarity is evidently 
satisfied at this level of accuracy. 
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C. Fundamental Scalar Interaction 



For the past 40 years, the weak interaction has been described by an equal mix of vector 
and axial-vector interactions that maximizes parity violation. The theory is known colloqui- 
ally as the 'V — A' theory. Despite the ever increasing precision possible in weak-interaction 
experiments, no defect has been found in the V — A theory. Prior to the establishment of the 
V — A theory, other forms of fundamental couplings, notably scalar and tensor interactions, 
were popular. Today there is still interest in searching for scalar and tensor interactions, not 
because we expect them to contribute importantly, but rather because we wish to establish 
limits to their possible contribution. 

A general form of th e weak-interaction Hamiltonian was written down by Jackson, 
Treiman and Wyld jlGS^ . In examining superallowed Fermi transitions, we are only in- 
terested in scalar and vector couplings, for which that Hamiltonian becomes 

Hs+V = {ijpiJn){Cs'(l)e<he + Cs'(j)el5(he) 

+ {i^pl^^i^n){Cv4>el^l(he +Cv'(f)^-f^,-f5(he)- (22) 

If we assume that the Hamiltonian is invariant under time reversal, then all the cou- 
pling constants must be real. Those coupling constants carrying a prime represent parity- 
nonconserving interactions. If we further assume that parity violation is maximal, then 
CI = Ci- In this limit, the scalar and vector terms can be written 

Hs+V = {^pCs^n){^ei^+l5)<h.) 

+ {^pCvl^^Pn) (0,7m(1 + 75)0fJ . (23) 

A nonrelativistic reduction of the hadron matrix element for the scalar and the time part 
of the vector interaction shows that they both reduce simply to constants in leading order. 
However, under charge conjugation the matrix element (ippCsipn) changes sign relative to 
{ippCvli'^^n)- Thus we write ±6*5 in the ensuing formulae with the upper sign being for 
electron emission and the lower sign for positron emission. The lepton matrix elements are 
different in the two terms in Eq. (j23|l so the contribution to the shape-correction function 
from the scalar interaction will involve a different combination of electron and neutrino 
radial functions than that from the vector interaction. The final formula for C{Z, W) is 
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+ ( rriKike, K) + Mxike, K) f 



Mxike, ky) + rriKike, k^ 



keW 

{mK{ke,k,) +MK{ke,k,))}, (24) 

where Mxike, kp) and mK{ke, k^) are the reduced matrix elements given in Eq. ()A13j) . which 
incorporate the radial functions, F{r) and /(r), defined in Eq. ()A14|) . The reduced matrix 
elements Mxike, ki) and fnK{ke, k^) are the same as Mxike, k^) and mK{ke, k^) except that 
the radial functions, F{r) and /(r), are replaced by F{r) and /(r), where 

F(r) = H{r) {-G^^ Jfe.-i(p.r) + Jfe.(p.r)} 
+D{r) {G+- jfc„_i(p^r) - G++ jkAP^r)} 

J(r) = h{r) Jk^APut) + G'_+ jkAVvr)] 

+d{r) jk.-iiPur) - G++jk^ippr)} . (25) 

The functions H, D, h and d are linear combinations of the electron functions, fni^) 
and Qnir), as given in Eq. ()A15|) : and the angular momentum factors G±± are defined 
in Eq. (IXT6|l . 



1. Order of magnitude estimates 

For a pure Fermi transition, the multipolarity of the transition operators is K = 0. 
Keeping only the lowest lepton partial waves, ke = I and A;,^ = 1, we expand the lepton 
radial functions in a power series in r. The order of magnitude of the lepton wave functions 
at small r are 

fi{r) = 1-... 
9~i{r) = I- ... 
/_i(r) = small 

gi{r) = small 
3o{Puf) = I- ... 

ji{Pur) = small, (26) 

We retain only /i(r), g-i{r) and jo{pur), setting their values to unity, and drop the other 
small terms. The angular momentum factors for K = L = s = are 6*++ = G = 1, and 
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G+_ = G_+ = 0. Then the amphtudes become 

Mo(l, 1) = CyMp + . . . 

mo(l, 1) = small, 
Mo(l, 1) = TCsMp + ... 

mo(l, 1) = small, (27) 

where Mp is the Fermi matrix element. The shape-correction function is then 

C{Z,W) = \Mp\'(c'y + Cl±'^-^CsCy) 

~ \Mf\^CI. (1 + bFji/W + . . .) , (28) 

where it is assumed that Cs ^ Cy- The term in bp'yi/W is called the Fierz interference 
term, with bp = ±2fiiCs/Cv- This is the well-known expression given by Jackson, Treiman 
and Wyld jl65|. Here /ii is one of the beta-decay Coulomb functions, Eq. ()A9|) . and is of 
order unity, and 71 = (1 — (aZ)'^Y^'^. 



2. Determining a limit on CsjCy 

With the results of our data survey, we can now search for any evidence of a \jW 
term in the shape-correction function, and hence set a limit ow Cs- The test is based on 
the corrected Tt values being a constant for all superallowed transitions between isospin 
T = 1 analogue states. For optimum sensitivity, we do not use Eq. for C{Z,W) in 
the evaluation of the statistical rate function, /, because of the extreme nature of some of 
the approximations made in deriving that equation. Instead we use the exact numerically 
computed expression. Since this calculated value of / depends on the value of Cs we simply 
treat Cs as an adjustable parameter and seek a value that minimizes iii a least-squares 
fit to the expression jFt = constant. The result is 

Cs/Cv = -(0.00005 ± 0.00130). (29) 

The sign of Cs/Cy is determined from the fit, since the calculated / depends on the in- 
terference between vector and scalar interactions. The interpretation of the sign is a little 
more delicate. We define C5 to be the strength of the scalar interaction in electron-emission 
beta decay, and this is the value quoted in Eq. (|^. Since all the superallowed data involve 
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positron emitters there is a sign change mentioned earher due to charge conjugation that 
operationally is included in the computations. The corresponding Fierz interference con- 
stant, bp, is just —2 times this quantity^: bp = 0.0001 ± 0.0026. Had we not assumed that 
parity was violated maximally then the outcome would be 



CsCy + CgCy 



-(0.00005 ±0.00130). (30) 



\Cv\^ + \C'v\^ + \Cs\^ + \C's\^ 
This result shows a factor of 30 reduction in the central value compared to our previously 
pubhshed result ^ , with the standard deviation being essentially unchanged. This is by far 
the most stringent limit on Cs/Cy ever obtained from nuclear beta decay. 

D. Induced Scalar Interaction 

If we consider only the vector part of the weak interaction, for composite spin-1/2 nucleons 
the most general form of that interaction is written 16(| as 

with being the four-momentum transfer, = {pp — Pn)^- The values of the coupling 
constants (vector), /m (weak magnetic) and fs (induced scalar) are prescribed if the CVC 
hypothesis - that the weak vector current is just an isospin rotation of the electromagnetic 
vector current - is correct. In particular, since CVC implies that t he ve ctor current is 



divergenceless, we anticipate that fs = 0. An independent argument jl67[, that there be 
no second-class currents in the hadronic weak interaction, also requires fs to vanish. Our 
goal in this section is to use the data from superallowed beta decay to set limits on the 
possible value of the induced scalar coupling constant, fs- This will provide a test of the 
CVC hypothesis and simultaneously set limits on the presence of second-class currents in 
the hadronic vector weak interaction. 



^ In our previous work described in ref. |J] and adopted in our subsequent publications, we explicitly included 
a minus sign in the formulae in recognition that all the superallowed Fermi transitions involved positron 
emitters. Thus the shape-correction function C{Z,W) was modified to C{Z,W){1 — ■jibp/W) and a fit 
of Tt{l — jibp/{W)) to a constant yielded a value of bp that was negative. Currently in Eq. we 
have defined bp such that C{Z, W) is modified to C{Z, W){1 + ''jibp/W) and hence we are now quoting 
bp with a positive sign. 
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1. Relation between fs and Cs 



Considering, then, just the induced scalar term in the vector part of the weak interaction, 

Hv{S) = ^,itfsq,)i^n 0e7/.(l + l5)(h., (32) 

we see that this term can be reorganised to match closely the Hamiltonian from the funda- 
mental scalar interaction shown in Eq. The momentum transfer, = {pp — pn)n = 
— {pe +Pue)^lJ can be moved into the lepton matrix element where, in combination with 7^, 
it can be replaced with the free-particle Dirac equation: 'y^{pe)^i4>e = irrieCpe, Jii{PuJi^i<Pue = 
imj;^(f)j:^, with me and mj7^ being the electron and neutrino masses, respectively. On setting 
the neutrino mass to zero, we find that Hv{S) becomes 

HviS) = ippTriefsi'n ^e(l + Ts)'/'^,. (33) 



This expression is equivalent to the fundamental scalar interaction in Eq. (j22j) with Cs 
simply replaced by rriefs- Thus, its effect on the shape-correction function can be des cribe d 
by the same replacement in Eq. ()28|). An equivalent result was obtained by Holstein 16S|. 



2. Determining a limit on fs 

We have now established the mathematical equivalence of the effects that fs and Cs have 
on the shape-correction function, C{Z,W). As a result, we can use Eq. (pUj) to conclude 
that 

mjs/gv = -(0.00005 ± 0.00130). (34) 

The sign of fs/gv follows the same convention as that described after Eq. ()29p. This result 
is a vindication for the CVC hypothesis, which predicts = ^ and /g = 0. Our result 
confirms this prediction at the level of 13 parts in 10^. As already mentioned, this result 
can also be interpreted as setting a limit on vector second-class currents in the semi-leptonic 
weak interaction, which therefore have not been observed here at the same level of precision. 



E. Right-hand Currents 

Let us no longer cons ider parity violation to be maximal. The general form of the weak 
interaction Hamiltonian jlGsl ] for just the vector couplings of relevance for superallowed beta 
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decay is 

Hv = {'ipplt^^n){Cv'(j)el^i<Pue + Cy^eT^Ts^F J (35) 

With Cy 7^ Cy we cannot associate the couphng constants with the hadron matrix elements 
as we did in Eq. ()23|) . Instead, the lepton and neutrino radial functions remain combined 
with Cy or Cy. The final formula for the shape- correct ion function then becomes 



C{Z, W) = E ^fce { ^ ( MliK, K) + ml{K, K) + NUK, K) + nlik,, k,) ) 

- ^ ( Mxike, ky)mK{ke, k^) + Nxike, k^)nK{ke, k^) ) |(36) 

where M^ike, k^), mxike, k^), NK{kf., k^), and nK{ke,k^) are reduced matrix elements as 
defined in Eq. ()A13|1 . with their respective radial functions being F{r), f{r), G{r), and g{r). 
These radial functions are 

F{r) = H{T){CyG^^j,,^^,M-C'yG^+j,SPur)} 

+D{r) {C'yG^^ Jk.^iiPur) - GyG++ 3kM} 
f{r) = h{r){GvG^^jk^^,{p.r)-Gi.G^^MMr)} 

+d{r) {GyG-+ 3K-i{,Pur) - GvG++ jK^Pur)} 
G(r) = i7(r){-C(.G__ j,_iM+CvG_+j,,.(p,r)} 

+D{r) jfc,_i(p.r) + G'yG++ jkAP^r)} 

g{r) = h{r) {-G'yG^^ JK-i{Pur) + GyG^+ Jfc.(p.r)} 

+d{r) {-GvG+^ 3K-~i{p.r) + JkM} (37) 

where the functions iJ, h and d are linear combinations of the electron functions, fnir) and 
gnif) as given in Eq. ()A15|) . The angular momentum factors G±^± are defined in Eq. ()A16jl . 



1. Order of magnitude estimates 

Consider a pure Fermi transition for which the multipolarity is = and only the lowest 
lepton partial waves, k^ = 1 and ky = 1, are kept. Then, as in Sect. IIV C ![ the amplitudes 
become 

Mo(l, 1) = GyMp + . . . 
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mo(l, 1) = small 

Noil, 1) = -C^Mf + ... 

rao(l, 1) = small, (38) 

The shape-correction function is then 

C{Z,W) = \MFm{c'y + C'y'). (39) 

We see that the dominant impact of the right-hand current is simply to scale the statistical 
rate function by ^(1 + Cy /Cy). This has no impact on the CVC test that demonstrates 
that J-'t = constant, but it does shift the value of the vector coupling constant and thus the 
deduced value of V^^. However, V^^ is obtained from the ratio of beta-decay to muon-decay 
rates, so before we can make any definitive statement on the effect of a right-hand current 
on V^j^, we must first examine the impact of that current on muon decay. We will show 
next that the correction due to a right-hand current is second order in small quantities in 
muon decay, but first order in vect or b eta decay. To this end we examine a more general 
Hamiltonian presented by Herczeg jlBflj ]. 



2. The effect on V^^ 

In the SU{2)l x f/(l) Standard Model, the semi-leptonic weak interaction Hamiltonian 
can be written schematically as 

HsM = ^VudiV -A){V- A), (40) 

where the first factor of V — A represents the lepton currents: V = (peln'Pv^ ^"^^ ^A = 
0g7^75077e, while the second V — A represents the hadron currents: V = ipplfii'n and —A = 
'(pp'JiM'Jb'ipn- The weak interaction coupling Gp/^/2 = g'^/8M^, where g is the basic coupling 
constant of the Weinberg-Salam Standard Model and is the mass of the exchanged 
ly-boson. 

, 17fll | considers an extension that is the most general form for non-derivative 



Herczeg 



local four-fermion couplings 



H = aLL{V -A){V-A) + aLR{V - A){V + A) 

+anL{V + A){V-A) + aRn{V + A){V + A), (41) 
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where again the first factor represents the lepton currents, the second the hadron currents. 
The lepton fields are now written as ^ = (pelfj.(Pv^ 'Pelt^'Pv^ depending whether the chirality 
of the neutrino is left-handed for V — A coupling or right-handed for V + A coupling. The 
neutrino states are in general linear combinations of mass eigenstates, 

= J:^UeiC|>f <I>1 = E^Ve^<P^ (42) 

where Uei and Vei are first-row elements of the neutrino mixing matrix. The observed beta 
decay probability is the sum of the probabiliti es of decays into the energetically allowed 



169L in assuming that the neutrinos 



neutrino mass eigenstates. We follow Herczeg 
produced in beta decay are light enough that the effects of their masses can be neglected. In 
particular, the terms that arise from the interference between amplitudes involving neutrinos 
of different chirality are dropped. Then the effect of neutrino mass mixing can be taken into 
account by our multiplying the coupling constants and a/,/? by ^/u^, and a^L and a^R 
by ^yv^ where 

Ue = E-\Ue^\' Ve = E'^\Ve^\'■ (43) 

The prime on the summation indicates that the sum extends only over the neutrinos that 
are light enough to be produced in beta decay. Note that if all the neutrinos are light for 
both left-handed and right-handed chiralities, then = Vg = 1 a.s a. consequence of the 
unitarity of the neutrino mixing matrix. 

Herczeg's Hamiltonian, Eq. ()4H1 . can be rewritten 

H = (aiL + aiR + ttRL + anR)VV 

+ {-0'LL - ClLR + ClRL + aRR)AV 

+ {-aLL + aLR - aRL + aRR)VA 

+ {aLL - aLR - aRL + aRR)AA (44) 



We can compare this with the Jackson, Treiman and Wyld (JTW) Hamiltonian jl65j . which 
in the current notation becomes 

HjTw = {CvV - C'yA)V + {-CaA + C'j,V)A 

= CvVV - C{,AV + C'^VA - CaAA (45) 
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Thus we identify the correspondences as^ 

Cv = ClLL + ClLR + CLrl + Orr 

C'y = O-LL + «Li?, — CLrl — ajiR 

Ca = —ClLL + ClLR + O-RL — O-RR 

C'a = -(^LL + ClLR - ClRL + ClRR (46) 

For Fermi beta decay, only the vector part of the weak hadron current contributes, so the 
decay rate, F^, as shown earher in Eq. (jH^ . is proportional to 



F^ OC i(|Cy|' + |C 



/ |2 

v\ 



I 1 2 I 1 2 

|«LL + 0,Lr\ + |«_RL + O-rr] 



(|1 + cilrI"^ + \aRL + aRR\ 
\aLL\Hl + 2ReaLR + ...) (47) 



where aij = aij/an 



To co ntinu e our determination of V^d we need to consider the purely leptonic muon decay. 
Herczeg jlGoj ] writes the effective Hamiltonian for muon decay in analogy to Eq. ()4ip as 

H = cll{V -A){V-A) + clr{V -A){V + A) 

+crl{V + A){V-A) + crr{V + A){V + A) (48) 

The coupling constants in Eqs. PHj) and (PT|) are related by the CKM matrix elements by 

oll = CllV^cI 

aLR = CLRe'^-y^d 
aRL = crlV^^ 

aRR = CRRe^'V^l (49) 

Here V^^ is the ud-mathx element of the CKM matrix for left-handed chirality quarks, 
and is for right-handed chirality quarks. The phase a is a CP-violating phase in the 



^ Herczeg |l69l ll7oj employs a metric that leads to a different sign on the 75 matrix, so his correspondences 



yield a different overall sign from ours for Cy and C'a. 
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right-handed CKM matrix. The decay rate, F^, for muon decay is proportional to 

OC \cll\^ + \clr\^ + \crl? + \crr\^ 

= \cll? (l + \CLR? + \crl? + \crr\^) (50) 

where = Cij/cLL- 

Combining Eqs. ()47p and ()50p. we obtain an expression that connects the ratio of beta- 
decay to muon-decay rates with the value of ^i^- 

l + \cLR?+\cRL\^+\cRRr ^ ^ 

In the Standard Model, only an and cll are non-zero; in any case, the quantities Uij and 
Cij with ij = LR, RL, or RR can certainly be considered small. The correction to the 
muon decay rate from right-handed interactions is therefore seen to be second order in 
small quantities, while the correction to Fermi beta decay rate is first order. Keeping only 
first-order small quantities, Eq. (jHTjl reduces to^ 

^ = \V,^f (I + 2ReaLR) . (52) 

If the neutrino masses are such that Me 7^ 1 and f e 7^ 1 then this equation is modified to 

where (we)/? in the numerator is given by the Ue in Eq. ()43p. with the sum extended over 
neutrinos light enough to be produced in beta decay, while in the denominator (ue)^ is 
given by the same expression but with the sum extended over neutrinos light enough to be 
produced in muon decay. Note that the Q-value for muon decay is 105 MeV, a factor of 
ten times larger than the Q-value for any Fermi beta decay we are considering. Also in 
Eq. is defined as J2[ \ Ui^i\^, where U^i are second-row elements of the neutrino mixing 
matrix. In what follows, we will assume (we)/? = (^^e)^ = i^^)^ = 1- 

Before proceeding to numeric limits, it is worth showing how the cur rent formulae relate 



to the simpler and more restrictive manifest left-right symmetric model 17ll |. In this model 



^ Herczeg \l7(\ also considers the possibility that the relation between purely leptonic and semi-leptonic 
Hamiltonians, Eq. 149|l . is not sufficiently general. He writes qll = {0'Ll)sm + cl'll^ with {aLL)sM = 
CLLV^d- Then Eq. ^ becomes T^/T^ = \V^J\l + 2Re{a'LL + aLR)), where a'^^ = a'LLl{aLL)sM- We 
will not pursue this further, but it is obvious the formulae above can accommodate this extension with a 
simple replacement of cllr with a^^ + cllr- 
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the departure from maximal parity violation is entirely due to the presence of a second 
ly-boson whose mass is much heavier than the usual W^-boson. If left-hand couplings are 
mediated by the boson, Wl, and right-hand couplings by Wr, then Wl and Wr will be 
linear combinations of the mass eigenstates Wi and W2, viz. 

Wl = W^i cos C + sin C 

Wr = e^'^(-VrisinC + 1^2CosC) (54) 

and u; is a CP violating phase. If it is further assumed that, apart from the different masses 
of the Wl and W2 bosons, the coupling constants and CKM matrix elements are identical 
for left-hand and right-hand couplings, then there are only two parameters in this model. 
These parameters are: 6 = (mi/m2)^ and (, where mi and 1712 are the masses of Wi and 
W2 respectively. Both parameters are small and, of course, are zero in the Standard Model. 
The parameters of Herczeg's Hamiltonian, Eq. ()41|). can be expressed in terms of 6 and ( 

— y -i/L — X 

All - - — 2 ^ud (^RR - ^ 

Smf 

aLR = aRL = -e^^C -C (55) 

for negligible CP- violating effects. In this limit, the expression for the ratio of Fermi beta 
to muon decay rates, Eq. (jH^ . reduces to 

^ = lKf.P(l-2C). (56) 

J- / 



This is the expression we used in our earlier work [l| to set limits on the extent of right-hand 
currents. 



3. Numeric Limit 



Let us now insert the experimental values from our survey data for the beta-decay and 
muon-decay rates to determine an experimental value for iKidP? which we will write as 
l^rflexpf This is the value we recorded earlier in Eq. (fT7|) . Then Eq. can be written as 

iKdIelpt = \V^,\\l + 2ReaLR) 

= (1 - \Vt? - K,?) (1 + '^ReaLR), (57) 
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where in the second hne we have inserted the condition for unitarity of the CKM matrix. 
Adopting the PDG's recommended values for \ V^g\ and (see text preceding Eq. (fT^) 
we obtain the following result from Eq. (j37j): 

(0.9482 ± 0.0008) = (0.9516 ± 0.0011) (1 + 2i?eaL/j) 

RettLR = -0.00176 ±0.00074. (58) 

Within the context of the manifest left-right symmetric model (see Eq. this result 
corresponds to ( = 0.00176 ± 0.00074, a similar value to the one we reported previously ^|. 
The result of a non-zero aiR or ( simply reflects the fact that the experimental values of the 
first-row CKM matrix elements do not satisfy the unitarity requirement. 

If, instead, we adopt the E865 value P] for V^s rather than the PDG average, then the 
result in Eq. (j58|) is modified to 

(0.9482 ± 0.0008) = (0.9484 ± 0.0014) (1 + 2i?eaL/j) 

RettLR = -0.00007 ±0.00084. (59) 

This result is consistent with no right-hand currents and unitarity being satisfied in the 
experimental CKM matrix elements. 



V. CONCLUSIONS 



Previous surveys of superallowed Fermi beta decay have at times noted disagreement 



among the derived J-'t values, and at other times agreement |^ |6[ . When disagreement was 
evident, subsequent attention paid to the problem led to both theoretical and experimental 
advances. As presented here, in Sect. IIIH the status now is of excellent agreement among 
all J-'t values - to better than 3 parts in 10^ over a wide range of nuclei from A = 10 to 
A = 74. Such agreement confirms the expectations of CVC, allows very restrictive limits to 
be set on the possible presence of scalar currents and makes it possible to go forward with 
confidence to the next steps - the determination of Vud and the unitarity test of the CKM 
matrix. 

The outstanding challenge at this time is that the value obtained for Vud, when combined 
with the current PDG-recommended values of V^s and Vub, leads to a unitarity test that 
fails by more than two standard deviations. There are no evident defects in the calculated 
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radiative and isospin-symmetry-breaking corrections that could remove this problem and, 
indeed, a shift in any one of these corrections large enough to restore unitarity would be 
almost impossible to justify [l|. Moreover, the derived value of Vud from nuclear decays 
has been remarkably stable for three decades despite a vast increase in the quantity of high 
quality data and many theoretical refinements in the calculations of the correction terms. 

So if any progress is to be made in firmly establishing (or eliminating) the discrepancy 
with unitarity, both theory and experiment must be brought to bear afresh on the principal 
sources of uncertainty. Although we will focus here on improving the nuclear contribution 
to the unitarity test, additional experiments are also required for neutron, pion and kaon 
decays. The first two provide independent, though so far much less precise, values for Vud] 
the third establishes the value of Vus, which may ultimately turn out to be solely responsible 
for restoring the CKM matrix to unitarity. Whatever the outcome for unitarity, however, 
the results of all these studies will provide crucial information, either in characterizing new 
physics beyond the standard model or in setting a tight limit on its possible existence. 

We have taken pains throughout this work to pay careful attention to all uncertainties, 
theoretical and experimental. In Section IIV Al we detailed the various contributions to the 
uncertainty in iKi^p. Of these, by far the largest is due to the nucleus-independent radiative 
correction, A^. Its uncertainty arises primarily from a box diagram involving the exchange 
of one W boson and one photon between the hadron and the lepton. To make the loop 
integration tractable, it is divided by a scale parameter in to hi gh- and low-energy portions. 
The high-energy contribution can be computed reliably but the low-energy one, as 

calculated originally by Sirlin |l73l |. depends on the choice of scale parameter. Sirlin chose 
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a reasonable range for this parameter, which has been retained by subsequent 



autho rs |174l . ll75| . It is this choice of range that drives the overall uncertainty on A^. Recent 
work |lj with effective field theories based on chiral perturbation theory has been unable 
to improve the situation: although this approach replaces the low-energy contributions to 
the loop diagrams by well-defined low-energy constants, the values of these constants are not 
known a priori. How to obtain a more refined, first -prin ciples computation of the low-energy 



172], but one of considerable importance 



contribution remains an open theoretical problem 
and urgency. Not only is this uncertainty the principal limitation on the precision with 
which Vud can be determined from nuclear superallowed (3 decay, but it will have a similar 
limiting effect on its determination from neutron and pion decays as well. 
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The next largest contributor to the error budget on \Vud\'^ is the isosp in-symmetry- 
breaking correction, 6c- Altho ugh uncertainties have been individually determined for the 

158l | of 6c — 6ns (see Table IIX|) . the dominant source of \Vud\'^ 



most recent calculations 



uncertainty attributable to 6c arises from the small systematic difference between the re- 



15S I , using 
159|, using 



suits from different theoretical techniques used to calculate 6c- Our approach 
Woods-Saxon functions, yields larger 6c values than the Ormand-Brown one 
Hartree-Fock functions. Here we have taken the democratic approach, considering that 
these two sets of calculations are equally likely to be correct and letting the difference be- 
tween their results determine a systematic uncertainty that we apply to the final result (see 
Eq. (Ha)). 

If reducing the uncertainty on must rank as the first priority for future theoretical 
work, then improving our confidence in 6c can be taken as the top priority challenge for 
experiment. Although there is no way to check the correctness of the absolute values of 6c 
from experiment, it is possible to check the nucleus-to- nucleus variations in the calculated 
values. The method, which is illustrated in Fig. |3J is based on the validity of the CVC 
hypothesis that the corrected J-'t values for the superallowed 0^ — » 0^ decays should be 
constant. In the figure we compare the uncorrected measured ft values (points and error 
bars) with the theoretical quantity jFt/ ({l + 6'fi){l — 6c + 6iys)) shown as a band, the width of 
which represents its estimated error. With the average J-'t value, jFt, taken from Table IIXI 
this comparison specifically tests the collective ability of all three calculated correction terms 
to reproduce the variations in ft from one transition to another. However, since 6'j^ is almost 
independent of Z when Z > 10, this test really probes directly the effectiveness of the 
calculated values of 6c — 6ns- 

It can be seen that there is remarkable agreement between theory and experiment. In 
assessing the significance of this agreement, it is important to recognize that the calculations 
of 6c and 6ns for Z < 26 are based on well-established shell-model wave functions and were 
further tuned to reproduce meas ured binding energies, charge radii and coefficients of the 
isobaric multiplet mass equation |l58| . The origins of the calculated correction terms for 
all cases are completely independent of the superallowed decay data. Thus, the agreement 
in the figure between the measured superallowed data points and the theoretical band is 
already a powerful validation of the calculated corrections used in determining that band. 
The validation becomes even more convincing when we consider that it would require a 
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FIG. 4: Experimental ft values plotted as a function of the charge on the daughter nucleus, Z. The 
bands represent the theoretical quantity J^t/{(1 + 6'f(){l — 6c + Sns))- The two groups distinguish 
those beta emitters whose parent nuclei have isospin = — 1 (darker shading) from those with 
= (lighter shading). 

pathological fault indeed in the theory to allow the observed nucleus-to-nucleus variations 
in 6c to be reproduced in such detail while failing to obtain the absolute values to comparable 
precision. Pleasing as the agreement in Fig. 0] is, though, new experiments can still improve 
the test, making it even more demanding, and can ultimately reduce the uncertainty in 6c 
further. 

These new experiments can follow different paths. One is to improve the precision on the 
nine superallowed transitions whose ft values are already known to within 0.15% or better. 
On the one hand, these are the easiest cases to study, all having stable daughters and all, 
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except ^"C, decaying predominantly (> 99%) through a single superallowed transition. On 
the other hand, they have been the subject of intense scrutiny for at least the past four 
decades and, given the number of careful measurements already published, the prospects 
for really significant improvements in these cases, at least in the near future, do not seem 
bright. Nevertheless, a glance at Fig. |21 shows that some modest improvements are certainly 
possible. If we accept as a goal that experiment should be a factor of two more precise than 
theory, then we see that the QEc-^alnes for ^°C, ^''O, ^^AP, ^^Sc and ''^V, the half-lives of 
ioq^ 2^A1™, ^^Sc and ^°Mn, and the branching ratio for ^''C can all bear improvement. 

A second path is to expand the number of precisely measured superallowed emitters to 
include cases for which the calculated nuclear-structure-dependent corrections are larger, or 
show larger variations from nuclide to nuclide, than the values applied to the nine currently 
best-known transitions. We argue that if the calculations reproduce the experimentally 
observed variations where they are large, then that must surely verify their reliability for 
the original nine transitions whose corrections are considerably smaller. The recent results 
for ^^Mg, ^^Ar and ^^Rb are the first cases of this type to reach sufficient precision that they 
can contribute to the test (see Fig. 0]) but more are sure to follow. We have included in our 
survey all cases that we believe are potential candidates within the next few years. 

Without doubt, these new cases present serious experimental challenges. In general, the 
parent nuclei are more exotic and thus more difficult to produce in pure and statistically 
useful quantities. They also exhibit more complex branching patterns, which for the cases 
with A > 62 include Gamow- Teller transitions that may be unobservable individually but 
collectively can play a non-negligible role 0]- These heavier nuclei also have very short 
half-lives, which currently limit the precision with which their Q^;^;- values can be measured. 
Even so, all these obstacles are obviously now being overcome and we may reasonably hope 
that before long there will not only be more cases with precisely measured parameters, but 
there will be more than one measurement of each parameter, an essential prerequisite for 
reliable results at the level of precision needed to constrain the correction terms. 

Although it would not impact significantly on the unitarity question, there is an additional 
reason to improve the precision with which the ft values are known, particularly for the cases 
with A < 26. A scalar current, if it exists, would manifest itself as a l/VT- dependence in 
the shape-correction function used in the /-value calculations (see Eq. (I2H1))- Since the 
superallowed transition energies decrease with A, this effect would be strongest for the 
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lightest nuclei, and ^ O. Improved precision in the !Ft values for those two nuclei would 
act directly to reduce the limits set on a possible scalar current. 

In conclusion, we can assert that world data for superallowed 0+ ^ 0"^ /3 decays strongly 
support the CVC expectation of an unrenormalizcd vector coupling constant, and also set 
a new limit on the possible existence of scalar currents. The nuclear-structure-dependent 
corrections used in the analyses of these data have so far stood up very favorably to ex- 
perimental tests, and the value currently obtained with them for Vud is deemed to be very 
robust, even though it is an important component of the top- row test of CKM unitarity that 
fails by more than two standard deviations. We have indicated the improvements required 
from both theory and experiment to increase the precision in future so as to produce a more 
definitive result for the unitarity test. 
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APPENDIX A: STATISTICAL RATE FUNCTION 

The statistical rate function is an integral over phase space, 

f = J^ pW{Wo-WfF{Z,W)S{Z,W)dW, (Al) 

where W is the electron total energy in electron rest-mass units, Wq is the maximum value 
of W, p = iW"^ — 1)^/^ is the electron momentum, Z is the charge number of the daughter 
nucleus (positive for electron emission, negative for positron emission), W) is the Fermi 
function, and S{Z, W) is the shape-correction function. If the shape-correction function is 
put to unity, the integral becomes the customarily defined one for beta decay, which we will 
denote as f^tsA- 

/stat = j"^" pW{Wo - WfF{Z, W)dW. (A2) 

The exact evaluation of / differs from /stat by 0.2% at A = 10 up to 5.7% at A ^ 74. Thus, 
to maintain 0.1% accuracy for / over that range, we must determine the shape-correction 
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function itself with 2% accuracy. Obtaining this accuracy requires consideration of the 
following issues: 

• The electron wave functions can no longer be simply those of the lowest partial wave 
(j = 1/2) generated by a point nuclear charge and evaluated at radius R, the nu- 
clear surface, but must be the exact functions for some chosen nuclear charge-density 
distribution; 

• The lepton wavefunctions exhibit some dependence over the nuclear volume, leading 
to what are called second-forbidden corrections. Furthermore, a more accurate treat- 
ment of the weak interaction leads to relativistic and induced- current corrections. All 
these effects must be incorporated since they impact on the nuclear matrix elements 
and inject a mild nuclear-structure dependence into the evaluation of /. 

• The atomic electrons cannot be ignored, but must be accommodated approximately 
in a screening correction; 

In what follows we describe the ingredients of a code we have written that incorporates 
these effects. It is based on the formalism of Behrens and Biihring jl66l| . Note that they 
define 

F{Z, W)S{Z, W) = FoLoC{Z, W), (A3) 

where Fq = 2F{Z, W)/{1 + 71). The purpose of their redefining the shape correction factor 
in this way was to remove the historic requirement to evaluate the electron wave functions 
at the nuclear surface. The product FqLq is given entirely in terms of the amplitudes of 
the electron wave function at the origin (see Eq. ()A9|) below), and C{Z,W) is the shape- 
correction function defined with respect to this choice. 

1. Electron Radial Wave Functions 

The wave function for the electron emitted in beta decay is given by the solution to 
the Dirac equation with an external electromagnetic field present, but restricted to the 
special case where the vector potential vanishes identically and the scalar potential is static 
and spherically symmetric. We solve in spherical coordinates and introduce a partial wave 
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expansion such that the basis states are written: 

^,^h^M-)Ur)x'-A^ (A4) 

where f^, Qk are radial functions, and Xk the usual spin- angular-momentum wave func- 
tions describing the coupling of the orbital angular momentum I and spin i to give a total 
angular momentum j with 2;-component /x: 



Xk = i^^l k^s\jfJ')Ylrni{r)Xms- (A5) 



minis 



The eigenvalue k, is 



and has values 



K = -j{j + l) + l{l + l)-l (A6) 

K = -(/ + l) = -(j + l) if J = / + i 

«: = / = (j + i) if J = / - i. (A7) 

The radial functions are solutions to a pair of coupled equations 

^+^-^9.-{W + l-V{r))f. = 
ar r 

^ - ^^/« + {W-1- V{r))g^ = 0. (A8) 
ar r 

Here V{r) is a spherically symmetric static potential that represents the interaction of the 
electron with the charge distribution of the nucleus. 

Our task is to solve the pair of coupled radial equations, Eq. ()A8|) . in three regions, 
< r < Ri, Ri < r < R2, R2 < r < oo, matching the solutions at each region boundary. 
The first region is the one over which the nuclear charge density is non-zero. Here, we 
establish a power-series solution, regular at the origin, as the starting solution and integrate 
numerically to Ri. In the second region, we have a pure Coulomb potential, for which an 
analytic solution can be found in terms of confluent hypergeometric functions of complex 
argument. The asymptotic solution in the third region is expressed in terms of the desired 
outgoing waves and a phase shift. The unknowns in the calculation are the phase shift, A^, 
and the normalization of the interior solution, a^'- they are determined from the matching 
conditions. 
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In the derivation of expressions for beta-decay observables, certain combinations of am- 
phtudes and phase shifts characterising the electron wave functions appear. These combi- 
nations are called the beta-decay Coulomb functions, and are generally organized so that 
they are of order unity, with corrections of order (aZ)"^. Since we are only interested here in 
the shape-correction function, the only such functions we need are FqLq, and yU^, where 
k = \k\, which are defined as 

_ a^-, + al, 
^o-i^o — 



2p 

,.,2 , .,2 



2 I ^,2 



= — — : — 2 ' (^9) 

where 7^ = {k"^ — (a ZY}y ^. Although these functions have actually been tabulated by 
Behrens and Janecke |l77 |. we have not used the tables but, in the interests of precision, 
have explicitly computed the functions from the calculated values of a^- 



2. Shape-correction function 

Behrens and Biihring 16(| give the following expression for the shape-correction function: 
C{Z,W) = Yl K{Ml{k,,k,)+ml{k,,k,) 



-MK{k„k,)mK{ke,k,)Y (AlO) 



keW 

where the sums over fcg and fcj, are partial-wave expansions of the electron and neutrino 
wave functions with ke = ie + \ and k^ = jy + \, and j,^ being the electron and neutrino 
total angular momenta. The integer K represents the multipolarity of the transition oper- 
ators and is limited to the range |je — ju\ ^ K < + ju- The functions M^ike, k^) and 
mx{ke,ky) are given in terms of form factors jl66| and we evaluate these form factors in 
the "impulse approximation". In this approximation, the nucleus is treated as a collection 
of non-interacting nucleons so it is only necessary to consider the weak interaction as acting 
upon a single nucleon. All the many-body aspects of the nucleus can thus be handled in the 
standard shell-model way. Let O be a one-body operator, which can be written 

O = J2{a\0\/3)aiap, (All) 
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where aj^ is the creation operator for a nucleon in quantum state a, and ap is an annihilation 
operator destroying a nucleon in state (3. The matrix element of O for many-body states 
becomes 

{f\0\z) = J2{a\0\P){f\aiap\t); (A12) 

that is, the matrix element is a linear combination of single-particle matrix elements, 
{a\0\P). The coefficients in the expansion, (/|a^a/3|i), are called one-body density matrix 
elements (OBDME). We leave it, then, as the job of the shell model to provide the OBDMEs 
and only deal here with the single-particle matrix elements. The functions Mx(fce, k,^) and 
mK{ke, ky) are now given in terms of reduced nuclear matrix elements for a single-particle 
transition jp ja by 



MK{k,,ky) = ^Y^i-Y'-'^ijo. \\F{r)TKLs\\jp) 

^ -Ji Ls 

rriKik,, ky) = ^ ^(-)^-^(j, II f{r)TKLs \\ Jp), (A13) 

-Ji Ls 

where j is a short-hand notation for (2j + 1)^/^ and Ji is the spin of the decaying nucleus. 
The radial functions are 

F(r) = i/(r){G__ 

f(r) = h{r) {G^^ jk^^iip^r) - jk,XPiyr)} 

+dir) Jk,.-iip.r) - G++ jk^r)} , (A14) 

with 

Hir) = HhSr) + 9-kSr)) 
D{r) = \{9kXr)-f-kXr)) 
hir) = \UkXr)-9~kXr)) 

d{r) = \{9kXr) + f-kXr)). (A15) 

Here /^(r) and ^'^(r) are radial electron functions^, while the spherical Bessel functions 
represent radial neutrino wave functions. The functions 6*+-, 6*-+ and G are 



^ The internal normalisations at the origin are here set to unity. Recall that these normalisations, have 
been separated out into the beta-decay Coulomb functions, FqLq, (see Eq. (lAQfl l. 
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short-hand notations for GKLsike^k^), GKLs{ke, -ku), GKLs{-ke,kt,), and GKLs{-ke, -ku) 
respectively, where the functions GxisiKe, Ku) contain all the angular momentum factors for 
the leptons: 



' K s L 

sKjejJJu {leOluO\LO) \ 1 



> . 



(A16) 



7 i / 

Lastly, the operators Trls depend on the angle and spin coordinates and are defined as 



rxL.(r) = (^0 + Ao)i''YLM{i)5K,L if s = 
= (V + A).z^Y;,iM(r) if s=l. 



(A17) 



where Vq, Aq are the time parts of the vector and axial-vector hadronic c urren ts and V, 
A are the space parts. Further 'Yxlm{'^) is a vector spherical harmonic 178l |. which in 
Eq. ()A17|) forms a scalar product with vectors V and A. 



3. Hadronic matrix element 

For nucleons, the vector and axial-vector interactions are written 

V^j, = gvj^ - fM<y^uqiy + ifsqti 

Were we discussing the weak interaction between point-like spin-1/2 fermions, then we would 
set (y^v = fl'A = 1, and /m = /s = /t = /p = 0. However, in considering nucleons, we 
recognize that they are not point-like and furthermore they are influenced by the presence 
of the strong interaction. Thus Eq. ()A18|) presents the most general form of a vector and 
axial-vector interaction that is consistent with Lorentz invariance and excludes momentum 
operators beyond the first power. Here = {pf — Pi)^ is the momentum transfer. The 
coefficients, in principle, could be functions of but, because of the low four-momentum 
transfer in beta decay, this dependence can be neglected and the coefficients are referred 
to as coupling constants with individual titles: being vector; g^^, axial- vector; /m, weak 
magnetic; /s, induced scalar; /x, induced tensor and /p, induced pseudoscalar. 
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Our aim is to reduce the matrix element ra/(V^ + A^)-Ui, where Uf and Ui are Dirac spinors 
characterising nucleons of momentum pj and pi, to the nonrelativistic form involving Pauli 
two- component spinors, 

xl,{Vo + y + Ao + A)xm,, (A19) 

by simply multiplying out the Dirac matrices involved and keeping terms to first order in 
|p|/Mn, and dropping terms in |p^|/M^ and higher order. This multiplication yields 

Vo = g. + fsiWo-V{r)) (A20) 

V = ^[p + t(T X q\ + U(T X q- f,q_ (A21) 

^0 = -^a.p-UcT.q (A22) 

A = -g^o- + U{Wo - V{r))cT, (A23) 

where p = P/ + Pi and q = P/ — Pi- Note that the large terms of order unity occur in Vq 
and A. Each of the four coupling constants denoted by an / are small and of order l/M^ 
and, as a consequence, terms in f /M^ have been dropped. Eqs. ()A20|) to ()A23|) are the 
quantities needed in the operators T^^^(r) in Eq. ()A17|) . 

4. Reduced matrix elements 



All the beta-decay observables 166| can be expressed in terms of the functions Mxike, k^) 
and mK{ke, k^) defined in Eq. ()A13|) . Here we are only interested in the shape-correction 
function, Eq. ()A10|) . which is a particularly simple combination of these quantities. Before 
proceding to evaluate Mxike^ky) and mxike, k^), we note that the expressions for both 
differ only in the presence of F{r) in one case and /(r) in the other. For simplicity in what 
follows, we will only explicitly display formulae incorporating /(r); obviously an equivalent 
set can be written with F{r). 

The operators T^^^(r) are shown in Eq. ()A17|) . but it is tidier if we incorporate the phase, 
(— )^~^ (see Eq. ()A13|) ). into the operator. All the operators (some after rearrangement) 
can then be expressed as a product of spherical tensors, one in orbital space and one in spin 
space. So, generically the operators take the form 

T^Lsir) = mn^'-'TKMiAL ® S^), (A24) 
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where is a spherical tensor in orbital space of multipolarity L, and S5 a spherical tensor 
in spin space^. We have introduced notation for a composite spherical tensor obtained from 
the combination of two other spherical tensors: viz. 

(A25) 

MlMs 

The single-particle wave functions from the shell model can also be expressed as products 
of orbital and spin space functions: 

\jm) = ^ {lmi\ms\jm)Rni{r)i^YimX^)Xms, (A26) 

m;r?is 

where n is the principal quantum number designating the number of nodes in the radial 
function. Notice the presence of 2' with the spherical harmonics^. 

The first step in evaluating the reduced matrix element is to factorize it into orbital and 
spin reduced matrix elements: 

{ila\)3a II TkLS II {lp\)3p) = 

{-y<-^A^LS)K {la II fir)AL II h){'^ II II i), (A27) 



where 



I i ? 



(A28) 



L S K 



Our conventions on reduced matrix elements are those of Edmonds jlTsj l. Next, we need to 
define two spin matrix elements, denoted Sq and Si and given by 

50 = (i II 1 II i) = V2 6s,o 

51 = (1 II II i) = Ve 6s,i; (A29) 
and two two orbital matrix elements denoted Ll and L(ji)l(Q). The first is 

Ll = {L II /(r)2^lL(r) || h) 



^ Note, the upper case S in Eq. l)A24p referring to the muhipolarity of the spin operator is not the same as 

the lower case s in Eq. l|IT7)l . 
® If one-body density matrix elements (OBDME) are imported from a shell-model calculation into this 

beta-decay environment, then it is important that these OBDME be computed with similar phases 

included in the definition of single-particle wave functions. 
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V47r I 



(A30) 



where the last factor is the radial integral: 

POD 

{RM{r)\Rp) = / Ra{r)f{r)R^{rydr. 

J 



(A31) 



The second involves derivative operators and requires a little more care. The matrix element 
is 



= {la II f{r)i'TL{Yj Q) II Ip), 



(A32) 



where Q is either p = p/ + pi or q = p/ — pj. Thus we write Q as p/ip, with the upper sign 
appropriate for p and the lower sign for q. Now Q is also —i^Vf ± Vj), where the gradient 
operator acts on either the initial or final nuclear wave function but not on the integrand, 
/(r). Thus the interpretation is as follows: {4>f\fQ\4>i) — — i{±(0/|/| V0i) — (V0/|/|0i)}. 
The result for iv(ji)i(Q) is 



-^(ji)l(Q) = i 



iu+i0+J+i^_y-L 



J 



X 



^U{lplU;lp-lL) 



-{-y+^-^U{lpJlal;la + lL) 



la J + 1 





la J Ip — ^ 





iaih + iy^'{Ra\f\{i-'i)R, 



iaih)'/' { Ra\f\ i i + ) R, 



\ 



Zq, + 1 J li3 





h{ia + iy^' {{i-'f)RM\Rp 



+ {-y+^-^U{lpJl,l;L-lL)\^^'' ^ ^ j^''^^''^'^' ((i^ + ^)^a|/l^/3/f. 

(ASS) 

where the upper sign is used for Q = p and the lower sign for Q = q. The [/-coefficient is 
a recoupling of three angular momenta and is related to a 6j-symbol: 



U{abcd; ef) = {-] 



■i+b+c+d 



ef 



a b e 
d c f 



(AS4) 
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Finally, we are ready to write down the specific reduced matrix elements for all the 
different hadronic components of the weak interaction, Eqs. ()A20|) to ()A23|) : 

Time-like Vector Current 

{m)ja II /(r)K)Z^Fx II {lp^)jf3) = {g. + fs {Wo + If)) A^K0)K Lk So. (A35) 

Space-like Axial Current 

{{laDja II /(r)A.z^Y^i II = (-)^-^ {-g, + {Wo + f f )) A^li)k S,. (A36) 

In Eqs. (IA35|) and ()A36|) . for simplicity we have replaced the function, V{r) (see Eqs. ()A20|) 
and ()A23|) ) by the potential due to a uniform charge distribution for small r, with replaced 
by its expectation value. In our computations we actually included the function V{r) in the 
integrand of the appropriate radial integral. 
Space-like Vector Current 

{m)ja II /(r)V.^^Y^i II = 

(-)^^^ l^-^^A(^KO)K -^{Ll)if (P) - fsA(^KO)K L(^Li)K{<i) Sq 

-V2J2 U{11KL; IJ) + /m) A^j,)k ^(Li)j(q) ^ij • (A37) 

Time-like Axial Current 



((/4)j„ II f{r)Ao^''YK \\ {hk)jp) 

J { gA 

J 



E(-)'"'' J {-^^(J^)^ kKi)j{p) Si - f.A^jDK L^kdAcO 5i} . (A38) 



5. Numerical Results 

The key ingredient for the computation of exact electron wave functions in beta decay is 
the charge-density distribution of the daughter nucleus. There are various parameterizations 
available in the literature, of which the following are the most common: 

• Two-parameter Fermi distribution (2pF). This charge density distribution, 

P{r) = — r-^, (A39) 

1 + exp|(r — c)/a\ 

has two parameters, c and a, other than its normalization. 

51 



• Three-parameter Fermi distribution (3pF). This is an extension of the two-parameter 
model, which introduces a dimensionless "wine-bottle" parameter, w, that impacts on 
the small-r behaviour of the density distribution. The functional form is 

• Three-parameter Gaussian distribution (3pG). This is an alternative three-parameter 
model with a Gaussian rather than a Fermi distribution: 

' l + exp{(r2-c2)/a2}- ^ ^ 

• Harmonic-oscillator distribution (HO). In light p-shell nuclei, where only s and p 
orbitals are occupied, a density distribution can be constructed from the harmonic 
oscillator radial wavef unctions. Its form is 

p(r) = po(l + arV&^) exp{-r^ /b"^), (A42) 

where b is the harmonic oscillator length parameter, and a is related to the number of 
p-shell protons, a = {Z — 2)/3. However, in practise both b and a are treated as free 
parameters and adjusted to fit the elastic-electron scattering data. 

These model distributions typically contain two or three parameters. Where possible, 
the parameters are determined from experimental data on elastic electron scattering, since 
the measured electron-scattering form factors are just the Fourier transforms of the charge- 
density distributions. A compilation of ch arge- density distributions determined from elec- 
tron scattering is given by De Vries et al. |l79|. We have assessed these data and selected 
for each daughter nucleus what we believe to be the 'best' value of the rms radius, (r^)^/^, 
and its probable error. In cases where data are not available on the isotope of interest, we 
have examined the nearest isotope that is available and applied a modest isotope shift to its 
value of (r^)^/^. Our final selected values are listed in Table 1X1 We also list the percentage 
uncertainty in the exact value of / due solely to the uncertainty in (r^)^/^. Clearly, the 
uncertainty in the charge-density distribution is not a factor in the determination of / to 
0.1% accuracy. 

Before the final evaluation of the statistical rate function, /, there are two further cor- 
rections to consider: for screening and recoil. To accommodate these corrections and to 

52 



TABLE X: Charge-density distributions from elastic electron-scattering data 179|. The radius 
parameter, in some cases, has been adjusted to reproduce (r^)^/^. Parameters, c and a, are in fm 
units, parameter w is dimensionless. 



Daughter 












A 

A/" 


Nucleus 


fm 


TV T _ J „1 b 

ivioaei 


c 


a 


w 


/o 


T — — r 














lOg 


2.45(10) 


HO 


1.709= 


0.837'' 




0.001 


14N 


2.52(2) 


3pF 


2.572 


0.5052 


-0.180 


0.000 


ISp 


2.90(3) 


2pF 


2.574 


0.567 




0.001 


22Na 


2.95(5) 


2pF 


2.750 


0.549 




0.001 


26 Al 


3.03(2) 


2pF 


2.791 


0.569 




0.001 


30p 


3.18(3) 


3pF 


3.350 


0.582 


-0.173 


0.002 


34C1 


3.39(2) 


3pF 


3.479 


0.599 


-0.100 


0.001 


38k 


3.41(4) 


3pF 


3.738 


0.585 


-0.201 


0.004 


42 Sc 


3.50(5) 


3pF 


3.794 


0.586 


-0.161 


0.004 


= 0: 














26 Mg 


3.06(5) 


2pF 


3.049 


0.523 




0.002 


34 g 


3.29(1) 


3pG 


2.810 


2.191 


0.160 


0.001 


38 Ar 


3.36(5) 


2pF 


3.590 


0.507 




0.004 


42Ca 


3.48(3) 


3pF 


3.765 


0.586 


-0.161 


0.002 


46 Ti 


3.61(3) 


2pF 


3.711 


0.588 




0.003 


50Cr 


3.66(4) 


2pF 


3.868 


0.566 




0.004 


54Fe 


3.69(2) 


3pG 


3.541 


2.270 


0.403 


0.003 


62 Zn 


3.90(2) 


3pG 


3.570 


2.465 


0.342 


0.005 


66 Ge 


4.04(4) 


2pF 


4.398 


0.585 




0.011 


70Se 


4.07(5) 


2pF 


4.442 


0.585 




0.011 


74Kr 


4.10(5) 


2pF 


4.489 


0.585 




0.013 



"Percentage uncertainty in / due to the uncertainty in (r^)^/^. 
''See Eas. dXal to (IX42ll . 

'^This is parameter, b, of Eq. (|A42|) in fm units. 
'^This is the dimensionless parameter, a, of Eq. ljA42p . 



remove - as is customary - the leading matrix element from the definition of /, we rewrite 
/ as follows: 

f = ^R{Wo) pW{Wo - WYFoLoC{Z, W)Q{Z, W)dW. (A43) 

Comparison with Eqs. ()A1|) and ()A3|) reveals three new factors, Q{Z, W), R{Wq) and S,. We 
will deal with them in that order. 
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The calculation of the Fermi fu nctio n presented so far makes no allowance for the screen- 
ing of the atomic electrons. Rose |l8c| was the first to find a simple analytic prescription to 
obtain the Fermi function for a screened field from the Fermi function for the corresponding 
unscreened field. That prescription is to incorporate a correction factor into the integrand 
for /: viz 

where W = W - Vq, p = {W^ - 1)^/^ and Vq = N{Z)a'^Z^/^, with Z being the electronic 
charge of the parent atom and N{Z) being a weak function of Z, which varies from = 1.42 



at Z = 8 to = 1.56 at Z = 29 (see Matese and Johnson [18l|). Since the factor Q{Z, W) 
yields a correction to / of order ~ 0.2%, we only need Rose's screening correction to be 
accurate to within 50 % of its central value in order to assure us an accuracy in / of 0.1%. 



Matese and Johnson 



18ll | have tested the Rose formula by comparing it with numerical 



solutions of the Dirac equation for a self-consistent Hartree-Fock-Slater potential. They 
conclude that Q{Z, W) has an accuracy of four significant figures or better for all energies 
except the very lowest in positron emitters. Since we integrate Q{Z, W) over the whole beta 
spectrum, which actually de-emphasizes the lowest positron energies, we conclude that the 
Rose formula has far more than sufficient accuracy for our purpose. 

The second new factor in Eq. ()A43|) is R{Wq), which is the correction for recoil: it 
recognizes that the daughter nucleus is not at rest but has a small amount of recoiling 
kinetic energy. As a result, the leptons' maximum energy is actually slightly less than Wq. 
The recoil correction |l66l | is 

3Wn 

where Ma is the average of the initial and final nuclear masses. For use in eq. ()A45|) . M4 
must, like Wq, be expressed in electron rest-mass units. The resulting correction is very 
small, being of order 0.02% for the superallowed beta decays from A = 10 to A = 74. 

Lastly, for allowed transitions it is customary to remove the leading matrix element 
from the definition of /. Thus, we have introduced ^ in Eq. ()A43|1 . where ^ = l/|Mirp for 
superallowed Fermi transitions, Mp being the Fermi matrix element. For pure Gamow- Teller 
transitions ^ = I/IMgtI"^, with Mgt being the Gamow- Teller matrix element. 

In Table IXll we list both the values of /stat, Eq- (jA2jl . and the exact values of / for cases of 
interest in superallowed beta decay. The relevant Qec value is listed as well. For the exact 
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TABLE XI: Comparison of statistical rate functions, /stat) /approx and the exact value, /. 



Parent 


keV 




/stat 


/ 


^/stat'* 
% 


Tz = -1: 












lOn 


191)7.9 


2.766b 


2.29484 


o onnon 


— U.ZD 


U 


2831. U 


4.5400 


/I O 1 A '~7 

42.6147 


/I o 7 /I o cr 

42.7485 


A O 1 

—0.31 


18 "NT^ 

JNe 


64:02. U 


5.6575 


133.867 


1 O yl yl O yl 

lo4.4o4 


A A£^ 

—0.46 


Mg 


41z4.D 


7.0716 


/lie OO/^ 

415.826 


A 1 O A Ar\ 

418.440 


A £?0 

—0.62 


26Si 


4836.9 


8.4656 


1014.75 


1023.28 


-0.83 


30s 


5459.5 


9.6840 


1945.49 


1967.05 


-1.10 


34 Ar 


6062.8 


10.8647 


3366.22 


3414.21 


-1.41 


38Ca 


6614.2 


11.9437 


5247.54 


5338.46 


-1.70 


42Ti 


7000.9 


12.7004 


6904.47 


7042.83 


-1.96 


Tz = 0: 












26t71 a 1 


4232.5 


IT OOOO 

7.2828 


474.691 


A 'TO 1 IT/^ 

478.176 


A 'TO 

—0.73 


34/^1 


5491.8 


9.7472 


1 ni'71 oc 

1971.95 


1 AA^ OA 

1996.39 


1 OO 

— 1.22 




6044.4 


10.8286 


OO AO AT 

3248.45 


OOAO 1 A 

3298.11) 


1 cr 1 

— 1.51 


42c!„ 
DC 


6425.7 


11.5748 


ylOr\1 '71 

4391.71 


yl yl TA yl 1 
4470.41 


— 1.76 


46 a /■ 


7050.7 


1 O 'TOTO 

12.7979 


7044.04 


Tl AA A/^ 

7199.96 


O 1 T 

—2.17 


50a /r 


7632.5 


13.9364 


1 A /I C/^ O 

10456.3 


1 ATOI 

10731.6 


o cr T 

—2.57 


54Co 


8242.7 


15.1305 


15281.4 


15750.1 


-2.98 


62Ga 


9171.0 


16.9472 


25187.9 


26247.6 


-4.04 


66 As 


9550.0 


17.6889 


30146.4 


31613.7 


-4.64 


™Br 


9970.0 


18.5108 


36605.4 


38602.2 


-5.17 


74Rb 


10416.5 


19.3846 


44606.1 


47277.8 


-5.65 



«A/stat = 100* (/stat-/)// 
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calculations, we imported one-body density matrix elements, OBDME, from a shell-model 
code. For each case we performed several shell-model calculations for various sets of effe ctive 
interactions and model spaces. We used, in fact, the same wave functions that we used jlSS^ 
to compute the nuclear-structure corrections 5c and 5]^s- Thus our / calculations can be 
considered to be entirely consistent with the calculation of the nuclear-structure-dependent 
corrections. The / calculation, however, is not very sensitive to the shell-model inputs. In 
light nuclei, different shell-model OBDME gave changes in / at the 0.01% level, increasing to 
around 0.1% m. A = 74, our heaviest-mass case. Where we have more than one shell- model 
calculation for a given nucleus, we have averaged the / values for the entry in Table IXIl 
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